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Janoes  A.  Mullen  and  David  Middleton 
Cruft  Laboratory,  Harvard  University 
Cambridge,  Massachusetts 


Abstract 


Nongaussian  noise  occurs  sufficiently  often  that  a study  of  its  points 
of  similarity  to  and  difference  from  gaussian  noise  is  desirable  . This  re- 
port considers  an  important  sub-class  of  non-normal  statistics,  viz.,  nearly 
normal  noise.  A new  form  for  the  nearly  gaussian  probability  densities  has 
been  found;  this  is  di.scussed,  compared  to  the  conventional  Edgeworth  series, 
and  related  to  earlier  work.  Next  a survey  of  physical  noise  sources  is  made 
in  order  to  classify  the  statistical  nature  of  the  noise  that  each  produces.  The 
rectification  of  nearly  normal  noise  and  a c-w  signal  in  a half-wave  vth- 
law  detector  is  analyzed  and  the  correlation  functions  of  the  output  obtained. 
The  behavior  of  the  output  is  compared  to  the  corresponding  results  for 
gaussian  noise  of  the  same  input  intensity.  For  the  linear  and  square-law 
detectors,  the  cases  of  prime  interest,  detailed  figures  and  a number  ofmore 
tractable  formulas  are  given.  Finally,  the  problem  of  finite  averaging  is 
briefly  considered. 


1 

Introduction 

This  report  considers  the  rectification  of  a-m  signals  and  nongaussian 
noise.  The  corresponding  problem  in  which  the  noise  is  gaussian  is  particu- 
larly important  because  most  noise  sources  are  gaussian  ones  and  because 
the  normal  distribution  is  easy  to  handle  analytically.  Quite  properly,  early 
work  considered  gaussian  noise  nearly  exclusively.  Nonetheless,  certain 
types  of  noise  do  not  have  normal  distributions,  so  that  an  analysis  of  the 
non-normal  case  is  desirable,  both  to  apply  directly  to  problems  where  it  is 
applicable,  and  to  help  to  indicate  in  uncertain  cases  how  critical  the  ar, sump- 
tion of  gaussian  statistics  is.  We  note  that  because  of  the  presence  of  a non- 
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linear  element  in  the  receiver,  viz.,  the  rectifier,  the  output  moments 
depend  on  the  entire  input  distributions,  so  that  the  input  statistics  affect 
even  the  simplest  of  the  output  properties. 

The  noise  of  the  present  paper  is  nearly  gaussian  so  that  the  proba- 
bility distributions  can  be  expanded  asymptotically  in  terms  of  the  corre- 
sponding normal  probability  distribution  and  its  derivatives.  This  is  the 
well-known  Edgeworth  series  [ 1]  ; however,  the  form  that  we  have  found 
most  useful,  which  involves  derivatives  with  respect  to  the  second  moments 
of  the  distribution,  has  apparently  not  been  mentioned  before.  Pearson 
used  the  idea  for  a special  distribution  in  order  to  facilitate  deriving  the 
usual  Edgeworth  series  and  Crofton  [ 3]  has  obtained  similar  equalities 
between  different  kinds  of  derivatives  in  a separate  derivation  of  the 
Edgeworth  series.  The  application  of  our  results  depends  partly  on  the 
fact  that,  after  passing  through  the  tuned  stages  of  the  receiver,  the  noise 
is  narrow-band.  We  have  obtained  relations  connecting  the  distributions 
of  the  envelope  and  phase  of  the  noise  wave  to  those  of  the  instantaneous 
value,  which  also  appear  to  be  new. 

A survey  of  physical  sources  of  noise  has  been  made  to  determi.jc, 

1 

whenever  possible,  whether  the  noise  produced  is  gaussian  or  nongaussian 
and  to  suggest  noise  models  for  use  in  the  analysis  of  the  detector.  The 
correlation  function  at  the  output  of  an  idealized  a-m  receiver  with  a 
vth-law  rectifier  has  been  found  when  the  input  is  a nearly  normal  noise 
and  a c-w  signal.  The  result  is  fairly  intractable  in  general,  but  formulas 
have  been  obtained.  Graphs  have  been  drawn  for  the  two  most  important 
cases,  linear  and  square -law  recti.'ier s , v/ith  various  representative  types 
of  noise  inputs.  The  general  effect  of  nongaussian  noise  compared  to 
gaussian  ■'ise  of  the  same  input  power  has  been  found  for  both  the  output 
power  and  the  output  correlation  function.  Finally,  the  effect  of  averaging 
the  output  over  a finite  time  has  been  examined,  again  qualitatively,  in  general, 
and  V, ,Th  explicit  results  and  figures  for  the  quadratic  detector. 
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II 

Probability  Distributions  of  Nearly  Gaussian  Noise 


2.  1 General  Properties  of  Probability  Distributions 

A sufficient  characterization  of  systems  containing  noise  requires  a 
statistical  description  of  the  noise,  viz.,  the  joint  probability  densities  of 
the  noise  at  a set  of  different  values  of  the  time,  or  alternatively,  know- 
ledge of  the  ensemble  from  which  the  sample  functions  actually  observed  are 
drawn  [4,5,6].  The  noise  processes  with  which  this  report  is  concerned 
are  Markoffian  [ 7,8]  , which  means  that  knowledge  of  the  joint  distribution 
at  two  times  is  enough  to  describe  the  process  completely.  Thus  two  distri- 
butions will  be  needed  for  each  variable. 

W^{y,t)dy  = the  probability  that  y will  lie  in  the  interval 
y,  y and  dy  at  time  t and. 


'^2^^1’h'  joint  probability  that  y will 

lie  in  y^^,  y^^  + uy^  at  nme  tj^,  and  that  y will  lie  in  y^> 
y2  + 2.  time  t^* 


As  the  time  difference  between  the  two  times  of  observation  is  increased, 
the  dependence  of  y^  on  y^  is  lessened  until  with  sufficient  separation,  yj^  and 
y^  are  independent;  that  is 


(2.1) 


lim 
t — >oo 

The  characteristic  function  of  the  noise  is  the  Fourier  transform  of  the 
probability  distribution  [9,10]  . 


i^iyi  + i\y^ 


(2.2) 


This  function  is  introduced  because  problems  are  often  simpler  when 
phrased  in  terms  of  transforms  than  when  stated  directly  in  terms  of  the 
probability  densities. 

A further  important  concept  is  that  of  stationarity  [ll]  , which  means 
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that  the  distributions  are  invariant  under  a linear  shift  of  all  the  times 
of  observation.  thus  does  not  depend  on  t,  and  depends  only  on 
the  time  difference 

An  important  set  of  parameters  associated  with  a distribution  .are 
its  moments,  defined  by 


(t) 
m n 


(2.3) 


Because  of  its  particular  importance,  is  given  a special  symbol, 

R(t),  and  a name,  the  correlation  function.  It  serves  as  a simple  but  not 
entirely  reliable  test  of  independence  of  the  noise  values  at  two  times  sepa- 
rated by  an  interval  t.*  Also,  the  spectral  distribution  of  the  noise  power 
can  be  found  from  the  correlation  function  by  using  the  Wiener -Khintchine 
theorem  [ 12,13]  , 


R(t)cos  wt  dt. 


W(f)cos  wt  dt. 


(2.4) 


When  all  the  Kth-order  moments  exist,  the  characteristic  function 
possesses  an  expansion  in  terms  of  the  moments  [9>10]  , 

<e^P{%Yi+  iS2yz)>av, 


K 


M-  (l) 
*^mn  ' 


Z. X n! 

M=0 

m+n=M 


+ o([5f 


(2.  5) 


A different  set  of  parameters  of  the  distribution  are  the  semi-invariants, 
defined  from  the  expansion  of  the  logarithm  of  the  characteristic  function 
[9,10]. 


❖ It  is  in  fact  a "wide  sensve"  test  for  independence,  in  the  sense  of  Doob,  sec  2.3. 
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oo  X 

inF2(5i«52*'^)  \ {nn  • (2.6) 

K=1 

k+l=K 

By  comparing  the  two  power  series,  relations  between  the  moments  and 
semi- invariants  can  be  found. 

The  sources  of  noise  considered  in  the  succeeding  section  are 
Markoffian  and  stationary;  furthermore,  they  are  representative  of  Poisson 
ensembles,  a particular  class  of  distributions  of  sums  of  independent  random 
variables , 

2.  2 The  Poisson  ensemble 

The  Poisson  noise  ensemble  is  composed  of  sums  of  independent 

variables  with  common  distributions  and  uniformly  distributed  times  of 

occurrence,  t.,  viz., 

J 

K 

j=l 

Physically,  the  basic  random  variables  represent  pulses  produced 
randomly  by  a noise  source.  These  hypotheses  imply  that  the  number  of 
pulses  occurring,  K,  is  a random  variable  with  a Poisson  distribution  [ 13,14]  , 
so  that  K,  as  well  as  the  set  of  t 's,  is  an  ensemble  parameter.  The  vari- 
ability of  K distinguishes  the  Poisson  model  from  the  random  walk  problem 
where  K is  a fixed,  known  number  [ 2,15]  . This  difference  is  not  a necessary 
consequence  of  either  model.  However,  with  superposed  pulses,  the  natural 
condition  of  uniformly  distributed  occurrence  times  leads  at  once  to  a distri- 
bution of  K.  In  a random  walk,  on  the  other  hand,  a fixed  number  of  steps 
seems  more  natural.  In  this  paper,  .''random  walk"  and  "Poisson"  models 
will  always  refer  to  ensembles  with  a fixed  or  variable  K,  respectively. 

The  distributions  of  V can  now  be  derived  following  Middleton  [16]  . 

Because  of  the  independence  between  pulses  the  characteristic  function 
for  V when  K pulses  occur  is  the  Kth  power  of  the  characteristic  function  for 


(2.7) 
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The  desired  characteristic  function  is  the  above  one  averaged  over  the  values 
of  K. 


2V 


oo  ^ 


K 


K=0 


(2.8) 


Here  y is  the  average  pulse  density,  the  average  number  of  pulses  per 
second  times  the  average  duration  of  a pulse.' 

The  above  equations  exhibit  a consequence  of  the  different  ensembles  for 
Poisson  or  random  walk  models.  In  the  Poisson  case,  the  semi-invariants  of 
V equal  y times  the  moments  of  v;  in  the  random  walk  (Eq,  2.  7),  the  semi- 
invariants  of  V equal  y times  the  semi -invariants  of  v,  i.e.,  the  fixed  value 
of  K is  y. 

In  the  general  model,  individual  pulses  are  themselves  random  variables 
with  distributions  of  amplitude,  phase,  duration,  and  occurrence  time,  viz.  , 

v(t-t')  = a h(p[t-t'];  r)  cos  [w^(t-t')  + ij']  • (2.9) 


Here  a is  the  amplitude  of  the  pulse;  h,  its  shape  factor,  which  depends 
on  the  occurrence  time,  t',  and  the  duration,  r(the  time  variable  has  been 
normalized  by  dividing  by  the  mean  duration  F = 1^  );  and  4<  is  the  phase  of 
the  pulse.  Most  systems  in  which  noise  is  important  are  spectrally  narrow - 
band;  accordingly  v(t)  has  been  written  as  a slowly  varying  envelope  factor, 
h(pt;r)  times  a rapidly  varying  "carrier"  wave.  In  terms  of  these  random 
variables,  the  characteristic  function  for  a single  pulse  is 


da  ^dr  ^ dr  ^di);  ^dv|;^dpt dpt'^dpt  j X w( a j ,a 2 , r j , r 2 ,t|;  j , vl>2 .Pt  j , pt '^) 


exp|ia  j5jh(P[t2-t'^] ; r ^) + ia2'52^(P[t-t'2]  J 


(2.  10) 


A specific  noise  source  will  not  usually  produce  pulses  variable  in  all 
these  quantities,  so  that  the  problem  of  evaluating  the  above  integral  is  not 
so  formidable  as  it  may  appear. 


The  moments  of  a distribution  are  always  of  interest;  in  this  paper, 
they  r re  particularly  so  since  we  shall  be  concerned  with  nearly  gaussian 
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distributions  where  only  the  first  four  sets  of  moments  appear,  not  the 
general  form  of  the  exact  distribution. 

For  simplicity,  let  us  consider  a noise  source  that  produces  pulses  at 
random  times  but  with  specific  amplitude,  shape,  and  phase.  Then  the  semi- 
invariants  of  the  output  are 


/ m+n 

\ a- 


i™(x;r)h’^(x  + pt;r)  cos^(wx/p-4j)cos^(a)^x/(3- vjj  + (j^t)d3^^. 


-oo 


(2.  11) 


where  ptj^  has  been  set  equal  to  a nev/  variable,  x,  and  t^  has  been  changed  to 
t + t. 

As  is  justified  in  Appendix  I,  because  of  the  narrow-band  structure  the 
phase  may  be  assumed  to  be  distributed  uniformly.  Carrying  out  the  phase 
average  first  then  simplifies  the  form  of  ^ considerably  in  specific 

cases.  In  general,  one  has 


\ 

\ 


mn 


(t)  = 0 


mn 


(t)  = X 


m + n odd. 


nm 


(-t). 


(2.  12) 


which  last  follows  directly  from  Eq.  (2.7)  above,  quite  independent  of  the 
narrow-band  assumption.  The  second  and  fourth  order  semi -invariants 
are  then 


Xll(t)  =v-2<ya  / h(x;r)h(x  + pt;r)d}^^os  o^t,  = i|)r^(t)  cos  o t = ij>r(t); 

^ 20  “ ^ 02  ~ ^ 1 1'^^  “ ^ ■ (2.13) 


oo 


1 /_4 


h^(x;r)h^(x + pt;r)dx>^^  (2  + cos  2o^t) , 


'-CO 


3 / 4 


r 


X3j^(t)  =Xj^2("l)  J li  (x;r)h(x^pt;r)dx^^^cos  o)^t, 

^oo 


TR189 


-8- 


[_^40  ^04  (2.14) 

It  is  convenient  to  separate  the  low- and  high-frequency  variations 

in  the  fourth-order  semi -invariants  just  as,  in  X.  ^(t).  the  v)jr^(t)  is  distinct 

from  the  carrier  part.  Defining  a low-frequency  factor  .A^  (t),  so  that 

A (0)  = X.  (0),  we  have 
^ 'mn'  ' mn'  ' 

\ 22(t)  = A22(t)[  2/3  + 1/3  cos  2w^t]  . (2.  15) 

3l(t)  = A3^(t)  cos  w^t 

^40  ^A4Q- 


2.  3 Nearly  gaussian  distributions 

A nearly  normal  distribution  may  be  asymptotically  expandedin  terms  of  the 
limiting  gaussian  form  in  inverse  powereof  y.  To  get  the  asymptotic  series,  we 
write  ( 2.6  ) as 


F2(5i.52-.t) 


1 + 


40 


4! 


^31^‘^  r , ^22^^) 
“Tl  5i  .?2  2!  2! 


2 

2 


+ 

X 


^13^^^  c 1e3  . ^ 04  '6  4 . ^ , -2, 
3!  ■^1^2'^  TP  ^2  ^ 

exp|-|[3f  + 2r(t)5i5^  +^21  • 


(2.16) 


The  symbol  0_(  y ^)  means  that,  after  the  transform  is  taken,  the 
^ _ 2 

neglected  part  is  of  order  y All  the  semi -invariants  are  of  order  y; 

is  of  order  so  that  is  of  order  y Thus  the  M = 4terms 

_ 1 ^ 

are  of  order  y ; the  next  succeeding  terms,  M = 6 and  the  square  ofM  = 4 
terms,  will  be  of  order  y ^ (odd  order  semi-invariants,  when  not  equal  to 
zero,  introduce  fractional  powers)  [ 17]  . The  probability  density  corret- 
sponding  to  the  above  characteristic  function  is 


W2(V^.V2;t) 


+ 


k3l(t) 

"Ti 


8 ^ ^ 22^^^ 
8Vj^8V2  8 
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V?  + - ZY^V^r(t) 


T\~^ 


3V^8V^ 


r “41 


04 

avi 


+ 0(y‘^) 


exp 


i-  Z ^Vl  - r^(t^ 


2it  ijj 


.VTT 


r^(t) 


(2.  17) 


This  is  the  two-dimensional  form  of  the  well-known  Edgeworth 
series  [ l]  . 

Another  form  of  this  probability  distribution  is  possible,  in  which 
derivativres  are  taken  with  respect  to  the  second  moments  (for  the  derivation 
of  this  result,  see  Appendix  II). 

The  characteristic  function  is 


+ 


+ 2 


+ y‘^)[ 


(2.  18) 


X expj^--^  [v|j^5j  + 2il^  r^(t)_5j52  cos  w^t  ^ ^ 


‘1^1  - 4^2  ■ 4> 


The  probability  distribution  can  be  obtained  from  this  by  transforming 
the  exponential  alone. 


The  series  (2.18)  was  first  found  by  Pearson  [ 2]  for  the  first  order 
density  in  the  random  walk,  using  the  properties  of  Bessel  functions;  but 
apparently  it  has  not  previously  been  extended  to  any  general  class  of  dis- 
tributions. The  only  condition  needed  is  that  the  random  variables  be  iso- 
tropic. Even  if  this  is  not  so,  a similar  series  involving  also  derivatives 
with  respect  to  the  first  moments  may  be  found.  Hov/ever,  the  latter  is 
sufficiently  involved  so  that  little  advantage  is  to  be  expected  in  using  it 
rather  than  the  Edgeworth  series. 


For  moments  after  a nonlinear  operation,  the  comparative  simplicity 
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in  the  analysis  when  the  noise  possesses  a gaussian  distribution  is  often 
a strong  point  in  favor  of  the  series  of  parametric  derivatives.  Applying 
the  differential  operator  to  the  result  for  ( nonstationary)  gaussian  noise 
is  likely  to  be  involved,  but  is  certainly  straightforward.  Furthermore, 
where  the  narrow-band  structure  of  the  output  is  important,  this  charac- 
teristic function  possesses  the  very  considerable  advantage  of  having  the 
high-frequency  part  appearing  only  in  the  exponent,  while  the  Edgeworth 
form  has  high-frequency  terms  of  different  orders  scattered  about  in  the 
various  semi-invariants  as  well  as  in  the  exponent. 

2.  4 Envelope  and  phase  distributions 

In  narrow  band  problems,  it  is  sometimes  more  convenient  to  take 
the  slowly  varying  parts  of  the  noise  wave,  as  the  random  variables,  rather 
than  the  instantaneous  value.  The  distributions  obtained  are  more  com- 
plicated, but  in  return  this  technique  entirely  eliminates  high-frequency 
terms  in  the  probability  distribution.  A narrow-band  ensemble  can  be 
written  as 

V(t;a,p)  = R(t;a)cos  [w^t  - 0(t;p)]  , (2.  19) 

wh€ire  the  envelope  and  phase  are  random  variables  not,  in  general,  in- 
dependent (which  means  that  the  ensemble  parameters  a and  p^re  function- 
ally related).  In  terms  of  in-phase  and  quadrature  components, 

V(t;  y [ a,p]  ,6[a,p])  = X(t;-Y ) cos  ujt  + Y(t;6)  sinwt.  (2.20) 

X and  Y are  not  usually  independent  either.  It  is  frequently  easier  to  find 
the  distributions  of  X 2Uid  Y and  afterward,  by  transforming  to  polar  coordin- 
ates, find  the  distributions  of  envelope  and  phase. 

From  the  equation  above,  it  is  evident  that  knowledge  of  the  distribu- 
tions of  V is  equivalent  to  knowledge  of  those  of  X and  Y.  The  explicit 
connection  between  these  two  different  representations  of  the  noise  wave  has 
been  found  under  the  condition  that  the  phase  is  distributed  uniformly.* 

*A  general  four-dimensional  distribution  has  thirty-one  non- zero  fourth- 
order  semi-invariamts  while  the  corresponding  distribution  of  V has 
only  five.  Finding  that  a connection  between  the  two  distributions  is  highly 
desirable  in  the  interest  of  algebraic  simplicity  alone. 


TR189 


-11- 


In  case  only  the  marginal  distributions  of  envelope  and  phase  changes  are 
needed,  the  condition  on  the  phase  distribution  may  be  dispensed  with.  It 
should  be  pointed  out  that,  while  the  fact  that  V is  narrow-band  and  nearly 
gaussian  implies  that  it  is  isotropic,  because  of  the  high-frequency  variations, 
that  is  not  true  of  the  slowly  varying  components,  precisely  because  the  high- 
frequency  part  has  been  removed.  Furthermore,  while  the  requirement  that 
the  phase  be  uniformly  distributed  is  satisfied  for  a good  many  noise  sources 
of  interest,  an  almost  equal  class  will  not  fit  this  scheme,  so  that  the  restric- 
tion is  more  severe  than  desirable.  Fortunately,  for  the  detection  problem, 
and  indeed  for  all  cases  where  only  the  low-frequency  output  is  needed,  the 
isotropy  condition  can  be  eliminated  by  working  with  the  distributions  of  the 
envelope  alone. 

The  high-frequency  variation  of  the  characteristic  function  of  V contains 
always  cosw^t.  To  put  the  high-frequency  dependence  in  evidence  let  us  write 
the  characteristic  function  as  cosoj^t).  Then,  if  the  charac- 

teftstic  function  for  X and  Y is  transformed  into  polar  coordinates. 


2X,Y 


(^Ci-t) 


F^x  y^5i  sin«i^,52COS  ^2'52sinf62lt), 


the  connection  between  the  two  is 


(2.21) 


^2X,Y^^1^°®  ^1’^1®^”^1’^2^°®  = ^2V^5i’52*'^’ 


(2.  22) 

a result  derived  in  Appendix  III,  on  the  assumption  of  a uniform 
phase  distribution.  Without  any  restrictions  on  the  phase,  however,  the 
distribution  of  envelopes  and  phase  difference  may  be  found.  We  have 


po  m Zn 


= R^R^ 


d^ 

(27t) 


Z 5,  52  dgjdS^ 


X + ^2  ^2 


(2.23) 


X ^2V^^r^2' 
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where  0 = Q2-0j^.l^  - 
As  an  important  special  case,  one  finds 

CD 

W^(R)  = R j ^J^(^R)F,^(^)d^. 
o 


(2.24) 


The  above  relations  are  true  of  narrow-band  distributions  in  general; 
these  must  now  be  applied  to  the  nearly  gaussian  distributions  which  are  our 
primary  concern.  If  there  is  an  unmodulated  carrier  present  as  well  as  the 
nearly  gaussian  noise,  the  instantaneous  value  of  the  sum  of  the  two  may  be 
written  variously  as 


Z(t)  = cosoi^t  + V(t), 

= [ A^  + X(t)]  cos  a^t  + Y(t)sin 
= x(t)  cos  oj^t  + y(t)  sin  cj^t. 


(2.  25) 


Using  Eqs.  (2.1),  (2.  21),  (2.25)  together,  the  distribution  of  the  envelope 
is  found  to  be 


W^(R)  = 


2A  R 
o 


-R^+A^ 

o 


4A  R , , 

-fr  t^o  + ^ ) 


2 ij; 


K—J 


(2.  26) 


^(R^+A^)  8(A^+R^) 


+ 8 


R^+  A^ 
o 


R 


2 ijj 


A nearly  gaussian  distribution  is  shown  in  Fig.  1 for  various  values 
of  signal-to-noise-power  ratio  together  with  gaussian  noise  of  the  same 
mean  square  value.  Before  the  nearly  gaussian  distribution  can  be  plotted, 
a value  of  the  fourth-order  semi-invariant  must  be  assigned.  The  graph 


FIG.  I ENVELOPE  PROBABILITY  DISTRIBUTIONS 
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here  represents  a noise  source  producing  square  c-w  pulses  of  constant  am- 

2 

plitude,  serving  as  an  approximation  to  clutter  in  a radar  system.  _A-^q/3!  ib 
then  becomes  l/4y  , where  y is  the  density  parameter  previously  defined 
(Eq.  2.8).  The  nearly  gaussian  distribution  has  more  small  deviations  from 
the  mode  and  a larger  tail  than  does  the  gaussian. 

As  pointed  out  previously,  the  random  walk  model  is  not  equivalent  to  the 
Poisson  one.  For  this  noise  source,  the  deviation  of  the  random  walk  distribu- 
tion from  the  gaussian  are  equal  in  sign  but  opposite  in  magnitude  to  the  devia- 
tion of  the  Poisson  model  distribution  in  the  figure  (for  equal  densitie  s) . In  general. 


^40 


_Ao 

,2 

Poisson 


Random 

Walk 


when  the  Poisson  and  random  walk  densities  are  the  same, 


(2.27) 


The  phase  distribution  presents  a more  difficult  problem.  The  sum  of  a 
sine  wave  and  narrow-band  noise  will  not  have  a uniform  phase  distribution  since 
the  phase  of  the  carrier  takes  a preferred  value.  If,  however,  the  noise  by  itself 
has  a uniform  phase  distribution,  then  the  desired  distribution  may  be  found  by 
transforming  coordinates  from  X,Y  of  Eq.  (2.20)  to  x,y  of  Eq.  (2.25). 

The  joint  distribution  of  envelope  and  phase  is  thus 


W^(R,0)  = 
whence 


r 

A 21 

f r2  + - 2A  Rcos  oI 

1+ 
1 + 

0 

o o 

04.2^ 

L ] 

(2.28) 


W^(0)  = 


, AtO  9^ 

94.2 


e 

2tt 


-p  sin^O 


j j 2 


or  as  a Fourier  series, 
oo 


W 


_ m/2 

X 7 € p cos  m« 

l(9)=  ^ 


-p(i^  + 1)  m + 1;  - p) 

A 


m! 


m = o 


(2. 29) 


40  m(m+2)  „ , m.  , _ , , . . 

qf F(~2-+2;  m + 1;  -p)V 


3!  4T 

The  envelope  distribution,  Eq.  (2.26)  can  also  be  derived  from  Eq.  (2.28) 
by  integrating  over  the  phase,  but  its  previous  derivation  directly  from 
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Eq.  (2.24)  emphasizes  that  no  assumption  on  the  phase  distribution  is 


necessary . 

The  second-order  distribution  is  needed  for  calculating  correlation 
functions.  For  a nearly  normal  noise  alone,  from  Eq.  (2.  23)  we  have 


"^40  r 8^ 


2 8 


r2(Rj.R2.ei.e,;t)  = R,R,|i+-3^[^4^J+3^^[A3i(t)4- 


.■^22^^V  1 ^ 8^  1 , ,,/  -2?i  f ■ ^‘^’*o^l^2^°®  ^■'■^1^2 

+— ^5 [-T — y + + Q(v  ) exp^ , J 


4)  8r^  '^1  ^^2  j [ 2(4jj4<2-4<  r^  ) J 

4>i=4>2=4> 

(2.  30) 

where  0 = 0^  - 8j,  and  the  parameters  are  those  of  Eq.  (2.  18). 

If  a sine -wave  signal  is  also  present,  the  distribution  of  the  envelopes 
and  phases  of  the  sum  can  be  found  by  transforming  variables  in  Eq.  (2.  25) 
to  the  polar  coordinates  corresponding  to  Eq.  (2.20).  The  result  for  a 
signal  and  noise  is 


^2^^  1’®1’^2’®2’  ^7 


= a + 


“^40r  8^  . 8^,  2 8 


84>j  8iJ< 


,^22^^^1  8^  . , 1 I Of  ■ 


J(2tt) 


® ‘h^  A^R  j(4i2-43r^)cosej+A^R^4ij-i|jrJcose. 

TT-J 5— j + t-t ^ 


) 2(4ji4<2-4j  ) 


v}ii  = V^ 


(2.31) 
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III 

Physical  Noise  Sources 


3.1  Introduction 

A search  through  the  literature  has  shown  that  the  Poisson  noise 
model  is  applicable  to  most  actual  noise  sources,  although  the  range  of 
densities  characteristic  of  nearly  gaussian  noise  is  found  only  occasionally. 
The  five  principal  groups  of  noise  sources  (not  all  mic roscopicklly  distinct) , 
considered  below  in  the  same  order,  are  those  associated  primarily  with 
(1)  radar,  (2)  sonar,  (3)  atmospheric  disturbances,  (4)  extra-terrestrial 
generators,  and  (5)  the  quantization  of  charge. 

3 .2  Sources  of  noise  in  radar 

Undesired  echoes  in  radar  are  an  example  of  Poisson  noise.  Besides 
the  desired  target  echoes  the  radar  picks  up  signals  from  the  surface  of  the 
sea  or  ground,  from  storms,  and  from  "window"  - metal  strips  used  to  pro- 
vide spurious  target  indications.  All  these  are  lumped  under  the  general  title 
of  "clutter"  [18,19]. 

In  a pulsed  radar  (with  a fixed  pulse  repetition  frequency),  the  input 
signal  following  a transmitter  pulse  represents  the  echo  obtained  from  a 
range  proportional  to  the  time  elapsed  since  theinitiation  of  the  transmitter 
pulse.  At  the  start  of  the  next  pulse,  the  range  is  set  back  to  zero  and  then 
increases  again,  as  shown  in  Fig.  2a  and  b.  In  analyzing  clutter,  we  want 
the  return  signal  from  one  fixed  range,  R^,  which,  as  shown  in  Fig.  2c, 
is  a number,  viz.  , the  amplitude  of  the  echo.  Before  being  presented  to 
the  observer,  the  signal  from  a fixed  range  is  converted  to  "boxcars," 

Fig.  2d,  and  then  smoothed,  as  in  Fig.  2e,  to  give  the  final  form.  The 
display  is  presented  on  a cathode-ray  tube  where  the  persistence  of  the 
phosphor  serves  as  the  boxcar  generator.  Variations  in  the  clutter  return 
between  successive  pulses  are  quite  small,  so  that  the  smoothed  represen- 
tation is  an  accurate  one. 

The  return,  nominally  from  a single  range,  R^.  actually  comes  from 
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a small  volume  around  within  which  targets  cannot  be  resolved.  The 
ncn-zero  angular  widths  of  the  antenna  beam,  © and^J.+  mean  that  targets 
in  a solid  angle  0<i>will  be  illuminated  simultaneously.  With  a pulse  du- 
ration, <r,  targets  in  a range  interval  of  length  c 't/2,  where  c is  the  pro- 
pagation velocity,  produce  echoes  which  overlap  the  return  from  R ; thus 

2 ° 
the  resolving  volume  is  R^  © ^ c t/2.  For  typical  radars,  the  pulse  duration 

is  about  1 psec  corresponding  to  a range  interval  of  150  yards  and  the  beam 

widths,  0.5-5°,  so  that,  up  to  ranges  of  1-10  miles,  the  shape  of  the  resolving 

volume  is  narrow  and  deep,  while  at  greater  ranges  it  is  broad  and  shallow. 

The  clutter  retdrn  at  the  detector  is  represented  by,  as  in  Eq.  (2.9), 


V(t^)  = 


(3.  1) 


where  the  scatterers  are  independent  and  uniformly  distributed  through- 
out the  resolving  volume.  Here  A is  the  amplitude  of  the  echo  from  the 
kth  scatterer  in  the  resolving  volume  and^jj^  its  phase;  e is  the  pulse  shape, 
a constant  of  the  transmitter,  but  also  a function  of  the  time  of  arrival  of 
the  received  pulse;  T is  the  duration  of  the  pulse,  also  constant;  and  is 
the  i-f  frequency.  The  average  number  of  scatterers  in  the  resolving  volume 
is  the  density  parameter  of  the  clutter,  because  the  ordinary  definition  of 
-y  as  the  average  number  of  pulses  per  second  times  the  mean  duration  of  a 
pulse  must  be  divided  by  the  amount  of  time  that  the  radar  receives  echoes 
from  the  single  range  R^. 

If  the  obstacles  are  moving  relative  to  the  radar,  the  clutter  at  a 
later  time  will  contain  variations  because  of  a shift  in  the  epochs 
tj^.  The  effect  of  motion  can  be  important;  however,  the  problem  is  treated 
in  Appendix  IV,  since  the  analysis  is  moderately  involved  and  does  not  bear 
directly  on  the  noise  density,  our  pri.mary  concern  in  this  section. 

Echoes  from  "window" 

One  important  type  of  clutter  is  "window,”  echoes  from  metal  strips. 
Analysis  of  the  qcho  from  a strip  [ 20]  shows  that  the  amplitude  depends  on 


*0  and  are  equivalent  rectangular  beam  widths  in  their  respective  directions. 


WAVEFORMS  FOR  A SINGLE-POINT  RANGE  DISPLAY 
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the  orientation  with  respect  to  the  direction  of  propagation,  but  that  the  time 
over  which  the  correlation  function  differs  from  zero  is  much  smaller  than 
the  rotation  time  of  a strip,  so  that  the  amplitude  return  is  effectively  con- 
stant in  time  i.  e.  , aj^  = a2).  The  total  number  of  strips  in  a resolving  vol- 
ume is  large  enough  so  that  the  noise  is  gaussian  if  the  strips  are  physically 
independent;  however,  they  might  cohere  in  clumps  so  that  there  are  only  a 
small  number  of  effectively  independent  scatterers.  Unfortunately,  present 
measurements  of  distributions  [ 21,22]  contain  enough  experimental  uncer- 
tainty to  mask  any  nearly  gaussian  deviations. 

Echoes  from  precipitation 

Rain,  hail,  and  snow  also  may  give  radar  echoes,  particularly  at  higher 
microwave  frequencies  [23]  . This  particular  problem  may  be  formulated  in 
terms  of  reflections  of  a plane  wave  from  an  aggregate  of  spheres.  Trinks[24] 
has  shown  that  for  Rayleigh  scattering  interaction  is  negligible  if  neighboring 
spheres  are  two  or  three  diameters  apart.  Meteorological  observations  of 
rainfall  [25,26,27,28]  have  shown  this  requirement  to  be  well  satisfied. 
Humphreys  [29]  gives  the  weight  of  rain  water  in  the  air  per  cubic  meter  as 
a function  of  the  rate  of  rainfall,  while  the  NEL  group  [27]  and  Laws  and 
Parsons  [26]  give  distributions  of  drop  sizes  as  a fxinction  of  rate  of  rainfall. 
The  drop  distributions  are  apparently  not  single -valued  functions  of  rain  in- 
tensity, but  they  are  consistent  enough  to  make  an  approximate  calculation  of 
the  density  parameter,  y , possible. 

As  the  rate  of  rainfall  varies  from  0.25-100  mm/hr.  (from  drizzle  to 

3 

cloudburst),  the  weight  of  suspended  water  goes  from  0.1-5  g/m  ; and  the 
average  drop  radius,  from  0.05-0.15  cm.  The  drop  distribution  gives  the 
relative  amount  of  water  at  the  ground  contributed  by  drops  of  a given  radius. 
This  is  the  important  distribution,  since  what  is  needed  here  is  the  amount 
of  water  in  drops  of  a given  radius,  not  the  number  of  such  drops.  The  dis- 
tribution should  be  corrected  by  dividing  by  the  terminal  velocity  distribution 
of  the  drops;  this  will  give  the  drop  distribution  in  the  air.  Since  the  larger 
drops  travel  faster,  this  correction  reduces  the  average  drop  size,  and  so 
increases  the  density  parameter.  As  the  lower  figure  is  found  to  be  large 
enough  to  provide  gaussian  noise,  additional  refinements  were  felt  to  be  un- 
necessary . 


TR189 


-18- 


The  average  number  of  drops  per  cubic  meter  is  then  150-400,  varying 
with  the  rate  of  rainfall.  This  corresponds  to  an  average  distance  between 
drops  of  12-8  cm.  Since  the  distUnce  required  for  negligible  interaction 
between  scattering  for  different  drops  is  0.3-0. 9 cm,  it  seems  clear  that  any 
lack  of  independence  in  the  scattering  is  small  enough  to  be  ignored  This 
many  drops,  scattering  independently,  give  a normal  noise  return  from  each 
cubic  meter. 

The  same  conclusion  may  be  reached  with  more  careful  reasoning, 
by  dividing  a cubic  meter  into  cubic  cells  such  that  the  scattering  from  pairs 
of  drops  is  interdependent  only  if  the  drops  are  in  neighboring  cellc.  The 
cells  prove  to  be  numerous  enough  (several  million  to  the  cubic  meter)  that 
the  possibility  of  having  more  than  one  drop  per  cell  is  negligible.  Since  the 
drops  are  distributed  uniformly  and  independently,  the  conditions  for  a 
Poisson  distribution  are  satisfied.  The  number  of  interacting  pairs  also 
possesses  a Poisson  distribution.  A calculation  of  probabilities  shows  that 
for  all  rain  intensities,  99.9  per  cent  of  the  time,  less  than  1 per  cent  of  the 
return  will  be  from  interacting  drops.  These  "exact"  figures  should  not  be 
taken  too  literally,  but  they  do  illustrate  the  overwhelming  likelihood  of  in- 
dependence. 

One  of  the  salient  characteristics  of  rainfall  is  its  inhomogeneity,  both 
temporally  and  spatially.  With  rain  gauges  less  than  5 meters  apart, 

Rado  [ 28]  found  variations  of  from  5 to  20  per  cent  in  the  measured  rainfall, 
while  the  NELi  group  [ 27]  , with  gauges  200  meters  apart,  found  much  larger 
differences. 

Thus,  in  one  resolving  volume,  the  rain  is  almost  certain  to  be  non- 
uniform.  If  the  non- uniformity  represents  merely  a statistical  fluctuation, 
then  unquestionably  the  noise  is  gaussian.  If  this  is  not  so,  on  the  other  hand, 
the  volume  may  have  to  be  split  into  subvolumes,  each  of  which  has  different 
statistical  parameters.  In  this  case,  the  subvolumes  will  be  larger  than  a 
cubic  meter,  so  that  the  total  return  will  be  a superposition  of  gaussian  noise 
from  each  subvolume,  and  therefore  gaussian  itself.  The  conclusion  from 
the  above  is  that  clutter  from  meteorological  precipitation  is  effectively  al- 
ways gaussian. 
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Echoes  from  the  sea  surface 

Another  source  of  clutter  is  sea  return— echoes  from  t;he  sea  sur- 
face [30,31]  . Exactly  what  reflects  the  incident  pulse  is  an  unsettled 
question;  quite  possibly  the  mechanism  may  differ  for  different  grazing 
angles  of  the  incident  beam.  There  are  three  possibilities;  (i)  specular 
reflection  from  the  waves  themselves,  (’ii)  reflections  from  ripples  on  the 
surface  of  size  comparable  to  a wavelength,  or  (iii)  reflections  from  spray 
droplets. 

The  radar  cross  section'?'  of  the  sea  for  high  grazing  angles  shows 
large  values  that  apparently  are  due  to  specular  reflection.  The  data  for 
low  grazing  angles  do  not  show  this  effect,  so  that  it  is  not  sufficient  to  ex- 
plain all  the  results  [31,3  2]  . 

The  dependence  of  cross  section  on  polarization  favors  the  droplet 
theory.  For  fairly  calm  seas,  the  cross  section  for  horizontal  polarization 
is  smaller  than  that  for  vertical;  the  difference  becomes  smaller  for  rougher 
seas  and  at  moderate  roughness  becomes  negligible.  This  can  be  explained 
by  destructive  interference  between  the  incident  and  reflected  waves  above 
the  surface  at  the  height  of  the  droplets.  A rough  sea  destroys  the  inter- 
ference pattern,  making  the  cross  sections  for  the  two  different  polariza- 
tions equal.  The  shifting  of  interference  patterns  with  grazing  angle  suffices 
to  explain  the  observed  dependence  of  cross  section  on  grazing  angle  [33]  . 

On  the  other  hand,  measurements  of  the  dependence  of  cross  section  on 

-4 

wavelength  [32]  , indicate  a X behavior,  which  favors  ripples  as  against 
both  droplets  and  specular  reflection.  Also,  the  drops  needed  to  produce 
the  observed  magnitude  of  sea  echo  are  much  larger  than  those  found  in 
practice.  The  approach  to  a limit  of  the  cross  section  at  low  wave  heights 
(2-4  ft)  also  favors  ripples,  since  with  high  waves,  more  and  bigger  drops 
and  steeper  wave  fronts  are  present.  The  last  two  factors  would  favor  a 


*The  radar  cross  section  of  an  obstacle  is  defined  as  the  area  intercepting 
that  amount  of  power  which,  when  scattered  isotropically,  produces  an  echo 
equal  to  that  observed  from  the  target.  For  sea  return,  the  cross  section 
is  expressed  as  a dimensionless  ratio  by  referring  it  to  unit  area  of  the  sea 
surface;  it  may  be  greater  than  one  if  the  actual  scatterers  are  not  isotropic. 
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continued  increase  of  cross  section  with  wave  height,  if  drops  or  specular 
reflection  were  the  primary  mechanism. 

The  nature  of  the  scatterers  is  not  settled;  the  observed  distributions 
are  usually  gaussian  to  within  the  limits  of  error  of  the  measurements. 

This  cannot  be  accepted  as  conclusive,  however,  because  of  the  major  effect 
of  deviations  from  the  normal  law  appear  in  regions  of  large  amplitude, 
which  is  the  part  of  the  distribution  most  difficult  to  obtain  accurately.  A 
long  period  of  observation  is  necessary  to  accumulate  enough  high-ampli- 
tude returns  to  decide  whether  the  clutter  is  gaussian  or  nearly  gaussian. 
Because  of  the  finite  lifetime  of  the  scatterers,  whether  drops,  ripples  or 
waves,  it  is  entirely  possible  that  the  clutter  statistics  be  non- stationary 
over  the  times  necessary  to  measure  the  tail  6f  the  probability  distribution. 
In  fact,  the  only  sea  return  distribution  in  Kerr  [34]  that  includes  the  larger 
amplitude  region  is  distinctly  non-normal,  from  whatever  cause. 

The  measured  relative  cross  section  is  -30  to  -70  db,  i.e.  , as  though 
- 3 -7 

10  to  10  of  the  illuminated  area  were  isotropically  reflecting.  This 

is  so  low  a value  as  to  argue  against  specular  reflection.  In  a typical  re- 

4 5 

solving  area,  there  are  10  or  10  sq.  ft.  , and,  assuming  a ripple  reflector 
to  cover  several  square  feet,  we  see  that  the  noise  may  well  be  nearly 
gaussian;  no  more  definite  conclusion  can  validly  be  drawn. 

3.3  Sources  of  noise  in  sonar 

Sonar,  considered  as  an  echo-ranging  device,  is  very  much  like  radar 
in  its  basic  operation.  The  wavelength  and  resolving  volume  are  of  the  same 
order  of  magnitude  as  in  radar,  because  of  compensating  changes  in  the 
other  parameters.  Typical  values  for  a sonar  [3  5]  are:  transmitter  fre- 

quency 24  Kc/s,  pulse  length  0.2  sec,  bandwidth  6-10°,  and  the  velocity 
of  sound  in  the  sea,  1600  M/s.  Then  the  resolving  area  is  about  300  yards 
in  range,  varying  laterally  with  range,  and  the  wavelength  is  about  6 cm. 
Serious  additional  problems  arise,  however,  because  of  the  dispersive 
properties  of  sea  water. 

Reverberation 

Sonar  clutter  comes  from  reflections  from  irregular  inhomogeneities 
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of  the  ocean,  and  from  the  irregular  boun^iaries  at  the  surface  and  bottom; 
it  is  usually  called  reverberation  rather  than  clutter  [36,37,38]  . The  im- 
portance of  the  three  types,  "volume,"  "surface,"  and  "bottom"  rever- 
beration varies  with  circumstances  as  shown  in  Table  3.1. 

Volume  reverberation  is  caused  by  sound  scattered  by  air  bubbles 
suspended  solid  matter  (e.g.  , plankton,  fish),  and  small  thermal  inhomoge- 
neities. The  measured  relative  cross  section,  i.e.,  scattering  area  per 

- 8 "5  “1 

unit  volume  varies  from  10  ° to  10  yd  at  a medium  range,  800  yds. 

This  makes  the  total  scattering  area  in  one  resolving  volume  (c  lO^yd^) 
.01-10  yd.  All  of  the  above-mentioned  scatterers,  except  for  fish,  are 
much  smaller  than  a wavelength,  so  that  the  individual  cross  sections 
should  be  very  small,  except  for  those  few  scatterers  of  a size  resonant 
for  the  particular  transmitter  frequency.  Then  the  received  noise  is 
"dense"  enough  to  be  gaussian. 

Table  3.1 

Predominant  Type  of  Reverberation 


long  range 
O 1 mile) 
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At  ranges  less  than  500  yards,  reverberation  from  the  surface 
may  be  20  to  30  db  higher  than  that  from  the  volume.  Surface  rever- 
beration drops  off  at  long  ranges  faster  than  volume  reverberation;  in 
addition,  it  depends  strongly  on  the  sea  state. 

As  a function  of  wind  speed,  the  return  starts  to  increase  at  8 mph, 
when  whitecaps  start  to  form,  and  continues  to  rise  until  it  saturates  at 
20  mph  with  the  reverberation  intensity  35  db  higher  than  the  initial  level. 
This  wind-speed  dependence  favors  bubbles  or  ripples  as  the  cause  of 
scattering;  however,  again  definite  conclusions  cannot  be  reached.  An 
analysis  of  the  scattering  from  a dense  layer  of  bubbles  shows  that  their 
cross  sections  are  not  large  enough  to  account  for  the  experimental  results 
by  themselves.  Scattering  from  ripple  patches  can  provide  the  necessary 
return  alone;  however,  the  observed  reverberation  is  probably  a combination 
of  the  two  [39]  . Bubbles  provide  a gaussian  return;  the  return  from  the 
ripples  cannot  be  calculated  without  some  knowledge  of  the  ratio  of  the 
geometrical  cross  section  to  the  scattering  one.  An  estimate  of  y is 
attempted  here  to  point  out  the  lacunae  in  our  knowledge  that  preclude 
definite  results. 

The  relative  cross  section  for  surface  reverberation  varies  from 

10  ^ to  10  At  300  yd,  the  resolving  area  is  9000  yd^,  which  cor- 

-3  2 

responds  to  9 to  9 x 10  yd  of  perfect  isotropic  reflectors.  The 

simplest  assumption  for  the  actual  reflectors  is  that  they  are  "perfectly 

diffuse,"  so  that  Lambert's  law  applies.  At  300  yd,  the  angle  of  incidence 

with  the  surface  is  about  6 degrees,  with  both  incident  and  reflected  flux 

corrected  to  normal  incidence,  900-0.9  yd  is  the  area  of  the  scatterers. 

Since  the  reverberation  intensity  is  found  to  be  frequency-independent,  it 

follows  that  the  scatterers  must  be  several  times  as  large  as  a wavelength, 

2 

say  1 yd  for  24  Kc/s  sound.  There  are  then  900-0.9  scatterers  per  re- 
solving area,  a density  range  including  both  gaussian  and  nearly  gaussian 
noise . 

The  above  figures  must  be  considered  as  simply  illustrative.  A 
proper  estimate  requires  knowledge,  rather  than  assumptions,  of  the  back- 
scattering  from  a ripple  and  the  average  size  and  structure  of  ripples.  In 
addition,  some  knowledge  of  bubble -density  is  needed  in  order  to  estimate 
how  much  of  the  surface  scattering  is  caused  by  bubbles. 
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Bottom  reverberation  is  important  only  in  shallow  water;  however, 
there  it  is  likely  to  be  the  dominant  factor.  This  scattering  arises  from 
irregular  features  of  the  bottom,  comparable  to  surface  reverberation  from 
ripples;  again,  not  enough  is  known  about  the  back  scattering  from  the  types 
of  obstacles  likely  to  be  encountered  to  make  a density  estimate  possible. 

Noise  from  marine  life 

A type  of  noise  that  has  no  counterpart  in  radar  is  that  produced  by 
marine  life  [40,41]  . The  important  distinction  is  that  not  echoes,  but  ex- 
ternally generated  noise  is  present.  Most  of  these  noises  are  markedly 
diurnal  and  seasonal;  one  important  type,  generated  by  the  snapping  shrimp, 
is  both  widespread  geographically  and  continually  produced.  Results  for 
this  noise  are  available  in  sufficient  detail  to  permit  an  estimate  of  the  ex- 
pected density  [42,43]  . 

The  shrimp  that  produce  this  noise  are  of  several  species,  all  of  which 
are  small  (about  two  inches  long)  and  possess  a large  snapping  claw,  which 
is  directly  responsible  for  the  noise.  These  shrimp  are  distributed  in 
tropical  and  subtropical  watersat  depths  less  than  60  meters,  with  a rock  or 
coral  bottom.  Investigation  off  San  Diego  has  definitely  located  the  shrimp 
as  the  source  of  crackling  noise  found  along  the  California  coast.  Among 
other  areas  in  which  both  snapping  shrimp  and  loud  underwater  crackling 
noise  are  found  are  Beaufort,  N.  C.  , Cape  Hatteras,  the  Bahamas,  Puget 
Sound,  and  several  locations  in  the  Central  and  Southwest  Pacific. 

One  of  the  important  features  of  this  noise  is  that  its  spectrum  is 

very  nearly  flat  from  1 to  15  kc,  while  most  underwater  noise  varies 

according  to  a (1/f)  law.  Shrimp  noise  is  thus  particularly  prominent  at 

the  higher  frequencies.  Measurements  on  an  individual  snap  show  an 

-4  . 2 

average  excess  pressure  of  60  db  above  2x10  ^yne  /cm  at  a range  of 

1 meter.  Corresponding  figures  from  four  shrimp  beds  found  off  San  Diego 
are  43  db  and  16  m,  34  db  and  30  m,  35  db  and  44  m,  and  45  db  and  20  m. 
The  density  can  be  found  as  the  total  power  divided  by  the  power  per  snap, 
both  measured  at  the  same  range.  For  the  four  measurements  above  , 

■y  is  5,2,6,  and  12  respectively, values  characteristic  of  the  lower  range  of 
nearly  gaussian  density.  The  noise  has  been  described  as  "like  static 
crashes"  or  "coal  going  down  a metal  chute,"  which  qualitati^'^ely  tends  to 
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indicate  a fairly  low  noise  density. 

3.3  Atmospheric  noise:  lightning 

Lightning  also  causes  noise;  however,  each  stroke  produces  only  one 
pulse  so  that  the  noise  is  not  dense  enough  to  be  nearly  gaussian.  Study 
of  lightning  shows  that  initially  there  is  a moderately  slow  downward  stroke 
followed  by  a rapid  upward  discharge,  of  durations  roughly  10  and  .05  Mc/s 
respectively.  After  this,  there  is  a further  slow  change  to  equilibrium. 
Superimposed  on  all  these  field  changes  are  small,  more  rapid  variations 
caused  by  fluctuations  in  the  resistance  of  the  ionized  path;  furthermore, 
multiple  reflections  from  the  ionosphere  may  occur  (44,45).  One  result  of 
the  presence  of  several  different  components  is  that  the  received  wave  form 
depends  strongly  on  the  distance  of  the  receiver  from  the  lightning  stroke. 

Near  to  the  stroke,  the  induction  field  will  predominate  and  a slowly  varying 
aperiodic  wave  form  will  be  observed;  far  away,  the  radiation  field  produces 
a quasi-periodic  wave  form.  The  various  components  all  belong  to  the  same 
stroke  and  different  strokes  are  spaced  widely  enough  so  that  the  noise  is  of 
low  density. 

A different  type  of  atmospheric  noise  comes  from  snow  squalls  [46]  . 

This  is  also  an  electric  discharge,  but  the  strokes  are  much  shorter  and 
weaker  than  in  thunderstorms.  While  the  strokes  are  more  frequent,  the 
essentially  distinct  character  of  the  individual  received  waveforms  shows 
that  this  noise  is  also  of  the  low-density  type. 

Precipitation  static 

Airplanes  flying  through  precipitation  in  the  atmosphere  become 
frictionally  charged.  If  the  storm  is  at  all  severe,  a corona  discharge  occixrs 
on  the  pointed  parts  of  the  aircraft,  e.  g.,  the  wing  tips  or  radio  antenna. 

The  noise  thus  produced  limits  radio  reception  considerably  and  may  cause 
a total  loss  of  intelligible  functioning.  Because  of  the  close  connection  between 
the  noise  and  the  weather  in  which  it  is  produced,  it  is  known  as  precipitation 
noise  [47,48]  . 

The  most  intense  noise  is  produced  from  corona  on  the  antenna.  Since 
aircraft  surfaces  in  practical  use  charge  negatively,  the  problem  is  that  of 
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negative  corona  from  a wire  [49]  . Initially  there  is  a glow  along  the  wire, 
which  with  increasing  potential,  collects  into  an  active  spot  giving  pulses 
from  a recurring  breakdown.  The  discharge  is  initiated  by  an  electron 
avalanche  in  the  high  field-strength  region  near  the  wire.  The  positive  ions 
left  behind,  moving  slowly,  eventually  fall  into  the  negative  wire  and  quench 
the  discharge  by  momentarily  lowering  the  potential  difference  across  the 
gap.  After  a recovery  time,  the  field  strength  becomes  high  enough  to  cause 
a discharge  again.  In  air  the  pulse  itself  lasts  about  0.1-0. 5 psec  while  the 
recovery  time  is  of  the  order  of  1 Mc/s.  These  recurring  pulses  are  called 
Trichel  pulses  after  their  discoverer  [ 50]  . 

As  the  field  strength  at  the  wire  increases,  the  frequency  of  the  Trichel 
pulses  does  also,  but  before  the  successive  pulses  merge  into  a continuous 
discharge  the  first  spot  saturates  and  another  active  spot  forms.  Because 
of  the  regular  recurrence  of  the  Trichel  pulses,  the  individual  pulses  do  not 
fit  the  Poisson  noise  model;  however,  the  sequences  from  different  spots  will 
do  so,  and,  as  seen  above,  more  spots  will  form  with  increasing  field  strength 
or  discharge  current. 

From  the  results  on  airplanes  [ 51]  , a check  with  the  corona  investi- 
gations is  possible  by  calculating  the  threshold  field.  At  the  corona  thresh- 
old, the  field  on  the  belly  of  the  test  airplane  was  250-300  volt/cm  which 
corresponds  to  voltages  of  25-30,000  volts.  For  an  antenna  0.04  inch  in  dia- 
meter, this  means  a field  at  the  surface  of  250-300,000  volt/cm,  which  is 
the  range  of  values  found  experimentally  [49.52]  . The  density  factor  may 
be  determined  from 

Y = (total  current)/(current/spot). 

The  total  current  figures  here  are  taken  from  the  results  of  Gunn  [51]  , but 
data  for  the  current  per  spot  are  not  available  from  the  results  of  Miller 
and  Loeb  [49]  , since  they  measured  the  total  corona  current  rather  than  the 
current  per  spot.  This  makes  a quantitative  determination  of  y i^npossible. 

In  air,  currents  start  at  0.1-1  pa  at  the  threshold  (this  is  from  one  spot)  and 
increase  to  values  of  50-100  pa,  by  which  time  multiple  spots  have  formed. 
The  total  discharge  current  from  the  airplane  will  be  from  1 pa  to  2000  pa, 
depending  on  the  severity  of  the  storm.  While  a calculation  of  y ^ function 


TR189 


-26- 


of  field  strength  is  not  possible,  it  seems  clear  that  the  noise  density, 
varying  continuously  from  low  to  high,  will  be  in  the  nearly  gaussian  range 
for  storms  of  intermediate  intensity.  This  is  qualitatively  borne  out  by 
aural  observations  of  the  character  of  the  noise  [53]  . 

3.4  Celestial  noise 

A further  source  of  noise  is  the  electromagnetic  radiation  from 
celestial  objects,  particularly  the  sun  [54]  Galactic  and  extra-galactic 
radiation  is  present  over  most  of  the  sky.  Since  the  sources  are  so  distant, 
the  mechanism  of  generation  is  not  completely  certain;  however,  it  appears 
that  the  radiation  can  all  come  from  point  sources.  The  continuous  radiation 
of  known  point  sources  indicates  that  the  mechanism  is  probably  a thermal 
one,  since  an  equilibrium  process  appears  necessary.  A thermal  source, 
which  here  means  a star,  contains  so  many  elementary  noise  generators 
that  the  noise  is  certainly  gaussian. 

Solar  radiation 

The  relative  nearness  of  the  sun  provides  the  possibility  of  deter- 
mining much  more  about  its  radiation.  Several  different  types  of  radiation 
in  the  meter-to-centimeter  bands  have  been  discerned,  one  kind  always  pre- 
sent, and  several  only  occasionally  present.  The  ever-present  radiation [55, 56] 
is  spoken  of  as  produced  by  the  quiet  sun,  since  it  provides  a minimum  solar 
noise  level.  As  it  is  always  present,  it  represents  an  equilibrium  process 
and  thus  is  thermal  and  gaussian. 

The  occasional  radiation,  produced  by  the  "active  sun,"  is  of  three 
types;  isolated  bursts,  outbursts,  and  enhanced  radiation.  Isolated  bursts 
occurring  for  a few  seconds  are  sometimes  present  in  solar  radiation  with 
peak  magnitudes  30-50  db  above  the  level  of  the  quiet,  sun.  The  bursts  are 
unpolarized  and  have  a sharp  but  smooth  rise  and  a slow  decay.  Their  in- 
frequent occurrence  marks  them  as  low-density  noise  [5  7]  . 

"Outbursts"  are  the  most  spectacular  of  the  various  forms  of  solar 
radiation.  Intensity  recordings  have  gone  off-scale  after  70  db  increases  in 
power.  Outbursts  last  for  several  minutes  at  a time;  they  are  apparently 
connected  with  visible  solar  flares~at  least  the  larger  of  them  are — and  are 
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often  the  cause  of  complete  interruption  of  short-wave  reception  Out- 
bursts do  not  seem  to  be  separable  into  superposed  pulses  and  so  can 
not  be  considered  as  Poisson  noise  [58,59]  . 

Enhanced  solar  radiation 

The  enhanced  radiation  from  the  active  sun  can  be  of  considerable 
importance.  While  not  the  strongest  form  of  solar  radiation,  it  is  30  db 
higher  (in  power)  than  the  quiet  radiation  and  lasts  longer  than  the  other 
active  forms  [60,61].  The  distinguishing  characteristics  of  enhanced 
radiation  are  its  duration  (several  days  rather  than  minutes  or  seconds,)  its 
circular  polarization,  and  its  correlation  with  sunspots.  Enhanced  radiation 
is  particularly  important  for  our  purposes,  since  it  is  the  only  type  which  ex- 
perimental results  show  to  be  nearly  normaL  At  least  two  widely  different 
origins  have  been  proposed  to  account  for  this  noise.  One  model  assumes  that 
the  radiation  is  due  to  the  thermal  motion  of  the  electrons  in  the  fields  near  a 
sunspot  [ 56,6  2]  . This  is  incoherent  and  depends  on  the  mean  thermal  energy; 
extremely  high  temperatures  (about  10^^  degrees  Kelvin)  are  required  to 
account  for  the  observed  intensities.  The  alternative  proposal  is  that  coherent 
oscillations  of  the  plasma  may  exist  [6  3>69]  . This  model  avoids  the  temperature 
difficulty,  but  must  explain  how  the  coherence  can  be  maintained  in  the  presence 
of  the  large  random  thermal  velocities  [63-67,69]  in  the  solar  atmospheres 
and  how  the  energy  can  radiate  from  the  generating  regions  where  the  group 
velocity  is  zero  [6  3,66-68]  . Opinions  are  too  unsettled,  and  detailed  know- 
ledge of  the  solar  atmosphere  near  sunspots  too  vague  to  make  a detailed 
discussion  of  present  theories  worthwhile  here;  nevertheless,  the  experimental 
results  show  the  features  needed  for  Poisson  noise  [5  ,70]  . 

Records  of  moderately  dense  noise  where  individual  bursts  can  not  be 
resolved  give  correlation  times  comparable  to  the  duration  of  individual 
bursts  discernible  on  quiet  days  [70]  . Wild  [60]  has  made  a careful  study 
of  enhanced  noise  over  the  frequency  band  70-130  Mc/s  for  a period  of 
several  months.  From  his  values  of  average  power  received  and  the  ob- 
served distributions  of  bursts  large  enough  to  be  individually  measurable, 
one  can  estimate  the  density  of  the  noise.  Wild  gives  an  experimental  in- 
tensity distribution  for  the  bursts.  The  middle  of  the  curve  follows  an 
exponential  law  fairly  closely;  at  the  low  end  the  pulses  were  occasionally 
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lost  in  the  background,  and  so  the  observed  number  of  burstsis  too  small; 

and  at  the  high  end,  the  infrequent  appearance  of  large  amplitude  bursts 

makes  the  sampling  error  large.  The  average  peak  power  in  a burst  is 

- 20  2 

found  to  be  2.5  x 10  watts/m  -c/s).  The  average  spectral  width  of  a 

burst  was  4.5  Me/ s between  quarter  maximum  points.  Since  the  spectral 

shape  was  found  to  be  essentially  gaussian,  the  width  of  the  equivalent 

rectangular -filter  with  the  same  peak  value  is  3.2  Mc/s.  Then  the  aver- 

- 20  2 

age  total  power  in  a burst  is  8 x 10  watts/m  . 

In  order  to  find  y,  only  the  total  power  received  needs  further  to  be 

calculated.  For  purposes  of  observation.  Wild  separated  the  received 

wave  into  a continuum  and  bursts  of  amplitude  10  watts  / (m  -c/s).  With 

the  observed  intensity  distribution,  the  average  value  of  the  continuum  is 

0.33  of  the  total  power;  including  the  60  Mc/s  bandwidth,  and  using  the  ob- 

- 20  2 

served  continuum  .75x  10  watts/(m  -c/s),  the  total  power  is  found  to  be 
-20  2 

135x10  watts/m  , wlience y equals  17,  in  the  nearly  gaussian  range. 

3 . 5 Noise  arising  from  the  quantization  of  charge:  shot  noise 

The  various  kinds  of  noise  caused  by  the  discreteness  of  the  elemen- 
tary particles  fit  the  Poisson  model.  Of  these,  shot  noise  is  the  simplest 
to  analyze.  In  a vacuum  tube,  the  current  contains  fluctuations  because  of 
the  discrete  nature  of  the  charge  [71].  As  long  as  each  electron  is  emitted 
from  the  cathode  independently  of  the  others,  viz.  , in  the  temperature-lim- 
ited and  retarding  field  regions,  the  Poisson  noise  requirement  is  satisfied, 

and  a density  estimate  is  straightforward.  The  transit  time  in  an  ordinary 
-9 

tube  is  about  10  sec,  so  that 

Y =<l>r/e  = 5x  10^°  <I>.  (3.2) 

The  noise  is  broad-band,  so  that  the  leading  correction  term  will  be  of  the 

-1/2  -4 
order  of  y '■>  ^ current  of  1 ma,  this  is  about  10  ; and  the  correc- 

tion does  not  rise  to  10  ^ until  the  current  drops  to  0.2  pa,  or  less. 

In  the  space-charge  region  of  the  characteristic,  the  electrons  are 
not  emitted  independently  and  a different  approach  to  the  density  problem 
is  demanded.  The  problem  now  is  that  each  emitted  electron  modifies  the 
space  charge  so  that  succeeding  electrons  find  a larger  potential  barrier  [72] 
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Obviously,  the  change  in  the  space  charge  is  produced  in  a fairly  small 
volume  around  the  individual  electron.  In  an  ordinary  receiving  tube  the 
potential  minimum  is  about  3 volts;  the  potential  due  to  the  moving  electron  is 
e/4-irr^,  sothat4xl0  cm  away  from  the  trajectory  of  the  emitted  electron,  the 
modification  of  the  potential  minimum  is  les  s than  1 per  cent  of  its  prior  value . 
The  quasi -stationary  potential  suffices  because  the  electrons  are  moving  in 
small  enough  fields  so  that  the  retarded  potentials  are  not  necessary. 

The  cathode  surface  may  be  divided  into  square  areas  of  such  size  that 
only  1 per  cent  of  the  electrons  emitted  in  the  area  will  produce  appreciable 
effects  outside  it,  i.e.,  a fringe  of  width  r^  around  the  edge  s of  the  square  con- 
tains 1 per  cent  of  the  total  area  of  the  square . Then  the  various  areas  are  ef- 
fectively independent  and  a noise  density  may  be  estimated.  Note  that,  since 
the  number  of  independent  areas  is  fixed,  the  random-walk  model  is  applicable , 

not  the  Poisson  noise  model . The  linear  dimension  of  one  of  the  se  areas  is 

-4  2 7 

400  r^,  or  1.6x10  cm.  With  a typical  cathode  area  of  0 . 5 cm  , y is  2x  10  , 

safely  in  the  gaussian  noise  region.  This  estimate  is  at  the  anode  of  the  tube  . 

If  the  network  intervening  between  this  anode  and  the  output  has  a rise  time  longer 

than  the  transit  time  of  the  electron  in  the  tube,  the  noise  pulses  will  be  stretched 

and  Y increased. 

A consistency  check  with  the  previous  re  suit  is  instructive  . Certainly 

the  effective  density  when  all  electrons  are  independent  will  be  greater  than  that 

when  they  are  not.  Infact,YT>  • /y  is  2x10  . The 

' 'Poisson  ‘space  charge  ^ ' 

space-charge  region  may  then  be  expected  to  start  with  currents  of  about  0 . 5 ma. 

A second  estimate  of  the  current  needed  for  a space-charge  effect  may 
be  made  as  follows:  if,on  the  average,one  electron  is  present  in  each  elementary 
area,  there  is  sufficient  overlapping  of  electron  emissions  so  that  the  as- 
sumption of  independence  is  untenable.  Only  about  0.1  of  the  emitted  elec- 
trons have  enough  energy  to  pass  the  potential  minimum,  and  consequently 

- 1 3 2 

the  charge  carried  from  one  elementary  area  is  6 x 10  coulomb/cm  . 

-9 

For  the  transit  time  of  10  sec,  a 1 cm  plate  to  cathode  spacing,  and  the 
cathode  area  0.5  cm^,  as  before,  this  is  a tube  current  of  0.3  ma,  in  order 
of  magnitude  agreement  with  the  previous  estimate. 

Photo -multiplier  noise 

Photo-multiplier  tubes  are  also  sources  of  noise  [ 73-75]  . The  im- 
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portant  distinction  between  these  and  other  tube  sources  is  that  the 
several  stages  of  electronic  amplification  produce  shot  noise  with  ele- 
mentary pulses  containing  many  electrons  instead  of  single  ones,  as  in 
the  ordinary  shot  effect.  This  means  that  the  photo -multiplier  noise 
current  fluctuates  more  intensely  than  ordinary  tube  current  of  the  same 
mean  value  [73]  . 

Besides  the  shot  effect,  noise  can  be  caused  by  spontaneous  emission 

from  the  photo  cathode  or  any  of  the  dynodes.  This  current  is  initially  very 

- 14 

small  (perhaps  10  a.  ),  but  electron  multiplication  again  magnifies  the 
effect  to  a level  requiring  consideration.  Such  emission  is  called  a "dark" 
current  since  it  flows  regardless  of  the  illumination  on  the  photo  cathode. 

Photo-multiplier  tubes  are  used  mainly  for  two  different  purposes, 
to  measure  light  intensities,  or  to  count  scintillations  in  radioactivity 
measurements . 

In  scintillation  counters,  radiation  causes  a fluorescent  screen  to 
produce  glowing  spots  that  activate  the  photo-multiplier  [76,77]  . 

The  signal  pulses  then  have  a duration  of  the  order  of  the  time  constant 
of  the  fluorescent  screen,  a few  microseconds.  The  noise  pulses,  on  the 
other  hand,  have  durations  of  the  order  of  the  transit  time  through  the  tube, 
about  10  sec  [74]  . The  input  time  constant  of  the  counter  pulse  ampli- 
fier is  made  equal  to  the  phosphor  decay  constant,  which  is  short  enough 
to  resolve  signal  pulses  and  long  enough  to  discriminate  against  the  noise 
pulses.  The  noise  pulses  at  the  anode  may  then  be  considered  6-functions 

since  they  are  of  such  short  duration  compared  to  the  input  time  constant. 

4 

In  a 931A  tube,  about  5 x 10  "dark-current"  pulses  are  emitted  per 
second  [75]  ; this  gives  a y of  about  0.1.  Dark  curi'ent  noise  is  thus  of  the 
low-density  type  in  scintillation  counters. 

For  light  intensity  measurements,  the  input  time  constant  is  large, 

(of  the  order  of  seconds,  in  fact)  in  order  to  smooth  out  fluctuations.  In 

5 

this  case,  the  dark  current  is  certainly  gaussian,  since  y is  about  10  . 
However,  the  signal  current  is  so  much  larger  that  the  noise  may  be  taken 
as  entirely  shot  effect,  with  the  dark  current  completely  neglected.  The 
transit  time  is  so  short  compared  to  the  input  time  constant  that  the  spread 
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may  be  neglected.  The  current  pulses  are  then 

i(t)  = ae  exp  - t/^rj  , (3.3) 

where  a is  the  overall  multiplier  gain. 

Now,  we  observe  from  Eq.  (2.  11)  and  (3. 3)  that 


(3.4) 


In  particular 


2<r 


(3.5) 


Notice  that  the  fluctuation  is  increased  by  at  least  a factor  ^a^over  the 
ordinary  shot  effect,  and  that  fluctuations  in  gain  produce  even  larger  out- 
put current  fluctuations . 


Since  X.  ^ equals  the.  mean  current,  we  have 


Y=<lX/<a>e,  (3.6) 

which  shows  that  y is  smaller  by  the  gain  of  the  multiplier  than  for  shot 

noise  of  the  same  mean  value  and  time  constant.  The  gain  of  the  tube  is 

5 6 

10  or  10  , so  that  this  effect  is  quite  important,  particularly  for  the  cur- 
rent at  the  photo -multiplier  anode , since  then  <r  is  the  transit  time  of  the  tube. 
In  the  input  circuit  of  the  voltage  amplifier  following  the  photo*- multiplier 
tube,  the  time  constant  is  about  1 sec,  so  that  y is  large  enough  for  correc- 
tion terms  to  be  safely  neglected. 

Thermal  noise  in  resistors 

Thermal  noise  from  resistors  is  one  of  the  commonest  sources  of 
noise.  From  a macroscopic  point  of  view  the  noise  is  sufficiently  well 
characterized  as  gaussian  with  the  well-known  power  spectrum  W(f)  = 

4k  TR  [ 78,79]  . The  noise  is  caused  by  the  random  motion  of  the  con- 
duction electrons  of  the  resistor,  achieving  its  randomness  by  collisions 
with  atoms  of  the  lattice  structure,  which  is  not  completely  regular.  For 
an  order-of-magnitude  estimate  of  the  noise  density,  the  Drude  theory  will 
be  sufficient  [80,81]  . In  this  simple  theory,  the  conduction  electrons  are 
independent,  which  is  necessary  for  the  Poisson  noise  model  to  be  directly 
applicable.  Furthermore,  the  mean  free  path  is  independent  of  the  electron 
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velocity,  while  the  scattering  from  a given  collision  is  isotropic. 

Here  y equals  the  mean  time  between  collision  times  the  average 
number  of  current  pulses  passing  through  a cross  section  per  unit  time. 

This  last  is  one  third  the  average  number  of  conduction  electrons  per  cubic 
centimeter  multiplied  by  the  product  of  the  mean  velocity  and  cross-sectional 
area.  Using  the  Maxwell -Boltzmann  velocity  distribution  and  the  independence 
of  the  mean  free  path,y  is  found  equal  to  4/3tr  times  the  product  of  mean  free 
path,  cross  section,  and  number  of  conduction  electrons  per  cubic  centimeter. 

One  has 

y = ^^NA,  (3.7) 


where  (j  is  the  mean  free  path,  N the  number  of  conduction  electrons  per  cubic 
centimeter,  and  A the  cross-sectional  area.  The  mean  free  path  may  be  ex- 
pressed in  terms  of  the  interatomic  distance,  d;  also,  the  number  of  conduction 

- 3 

electrons  is  proportional  to  the  number  of  atoms,  which  should  vary  as  d 
We  have 


A ; 


(3.8) 


• 8 

d is  of  the  order  of  10  cm.  , k^,  about  100,  and  kj^  of  the  order  unity,  so  that 

y = 4x  10^^  A(cm^)  (3.9) 


The  leading  correction  term  to  the  normal  distribution  is  of  order 
y which  is  certainly  negligible  here. 


Non-equilibrium  sources  of  noiee 

Besides  thermal  noise,  conduction  in  solids  causes  current  noise, 
semiconductor  noise,  and  flicker  noise.  These  last  are  nonequilibrium 
effects;  their  spectrum  is  directly  proportional  to  the  mean-square  current 
and  inversely  proportional  to  frequency.  These  phenomena  are  complex, 
each  depending  apparently  on  several  different  mechanisms  whose  relative 
importance  is  critically  dependent  on  factors  not  ordinarily  vinder  control. 

Flicker  Noise 


Flicker  noise  is  caused  by  variations  in  the  emission  from  thermionic 
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cathodes  over  times  large  compared  to  a single  transit  time  [82,83,84] 

For  some  reason,  certain  regions  of  the  cathode  surface  emit  particularly 
strongly  for  short  times.  This  corresponds  to  fluctuations  in  the  work 
function,  probably  caused  by  diffusing  ions.  A diffusion  process  leads  to  a 
1/(1  +a?<r^  spectrum,  where«7"is  the  lifetime  of  a diffusing  ion;  however,  if 
an  appropriate  distribution  of  lifetimes  is  assumed,  a l/cj  spectrum  can  be 
obtained  over  any  finite  interval  [84]  . At  present  knowledge  of  the  mechanism 
is  insufficient  to  warrant  an  estimate  of  noise  density. 

Excess  noise  in  solids 

Certain  types  of  resistors  develop  noise  in  excess  of  the  Nyquist  for- 
mula when  current  is  passing  through  them  [83-85]  . Besides  the  dependence 
of  noise  power  on  the  construction  of  the  resistor  (e.g.  , whether  it  is  carbon 
composition  or  "metallized”),  different  resistors  of  the  same  kind  of  con- 
struction and  of  the  same  value  of  resistance,  can  give  noise  powers  differing 
by  a factor  20.  Here  again,  exactly  to  what  the  fluctuations  should  be  as- 
cribed is  not  known  in  enough  detail  to  allow  a calculation  of  y . 

In  semiconductors,  two  types  of  current  carriers  are  present,  each  of 
which  produces  noise.  Excess  noise,  with  the  characteristic  1/f  spectrum, 
is  apparently  caused  by  fluctuations  in  the  number  of  minority  carriers  [86]  . 
Minority  carriers  affect  the  conductivity  strongly,  since  they  add  to  the  total 
number  of  current  carriers  themselves  and  by  lessening  the  space  charge, 
also  increase  the  number  of  majority  carriers  [87]  . In  the  material  not  near 
an  electrode,  the  current  density  will  be  nearly  uniform;  using  the  same 
reasoning  as  for  thermal  noise,  we  see  that  y ^ ■£,  NA.  The  cross-sectional 
area  is  about  the  same  as  for  thermal  noise;  and  the  mean  free  path  is  com- 
parable [88]  ; however,  the  concentration  of  minority  carriers  is  much  less, 

10  or  10  of  the  total  number  of  atoms  [89]  . The  excess  noise  is  thus 
much  less  "dense"  than  thermal  noise;  nevertheless,  y will  be  sufficiently 
large  so  that  the  noise  is  effectively  gaussian.  Evidence  has  been  presented 
that  another  component,  due  to  majority  carriers,  is  present,  becoming  rela- 
tively more  important  at  the  higher  frequencies  (above  10  kc/s)  [90,91]  . 

Since  this  noise  is  denser  than  minority  carrier  noise,  it  will  be 
gaussian  also. 
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Barkhausen  noise 

The  noise  in  ferromagnetic  materials  caused  by  changes  in  the  domain 
structure  [92]  is  called  Barkhausen  noise.  As  a magnetic  field  is  applied, 
two  different  processes  take  place.  At  low  values  of  the  impressed  field,  the 
domain  boundaries  move,  with  the  result  that  favorably  oriented  domains  be- 
come larger  at  the  expense  of  slightly  misoriented  domains.  With  greater  im- 
pressed fields,  the  direction  of  magnetization  of  domains  originally  oriented 
nearly  perpendicular  to  the  field  is  rotated  toward  alignment  with  the  field, 
without  any  wall  motion.  Domain  rotation  produces  noise  because  of  the  dis- 
crete character  of  the  changes  of  magnetization  [93,94]  . Inhomogeneities 
making  the  wall  motion  irregular  may  also  produce  noise  [95]  . As  the  domain 
wall  moves  past  an  imperfection,  subsidiary  walls  reaching  from  the  irregu- 
larity to  the  main  front  are  formed;  with  continued  expansion  of  the  main  wall, 
the  subsidiary  walls  stretch  imtil  they  break,  giving  a smooth  main  wall  and 
leaving  a small  stationary  domain  around  the  imperfection.  Barkhausen 
noise  has  been  observed  when  no  domain  rotation  was  occurring,  favoring  the 
above  mechanism.  Quantitative  statements  about  the  relative  importance  of 
these  two  mechanisms  or  estimates  of  the  noise  density,  do  not  yet  appear 
pos  sible . 


FIG.  3 BLOCK  DIAGRAM  OF  AN  A-M  RECEIVER 
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IV 

Pvectification  of  Nearly  Gaussian  Noise 

The  next  task  is  to  use  the  distributions  previously  found  in  sectionll, 
with  specific  fornmsofthe  correlation  function  and  semi-invariants  suggested 
by  section  III,  to  determine  the  most  important  statistical  parameters  of  the 
output  of  an  a-m  receiver  with  a c-w  signal  and  nearly  gaussian  noise  im- 
pressed on  the  input. 


4.  1 Model  of  an  a-m  receiver 

In  an  amplitude  modulation  system,  the  signal  is  transmitted  from 
source  to  receiver  at  a relatively  high  frequency,  which  is  demodulated  by 
the  receiver  to  recover  the  envelope  of  the  high-frequency  variation,  which 
represents  the  original  signal.  The  noise  is  generated  in  the  channel  and 
in  the  first  sections  of  the  receiver.  The  receiver  contains  a band-pass 
filter  which  eliminates  all  spectral  components  except  those  in  a narrow- 
band  around  the  central  frequency,  and  a demodulator,  composed  of  a half- 
wave v-th  law  detector  followed  by  a low-pass  filter.  The  relation  between 
output  and  input  of  a nonlinear  device  is  called  the  dynamic  transfer  charac- 
teristic; for  the  half-wave  v th-law  device,  we  have 


I = g(V) 


pv''  V>  0 
0 V<  0 


(4.  1) 


where  I is  the  output,  V the  input,  and  p an  appropriate  proportionality 
constant.  The  input,  V,  is  a narrow-band  signal  and  narrow-band  noise , 
but  the  nonlinear  device  changes  the  spectral  composition,  so  that  the  out- 
put is  composed  of  spectral  zones  lying  at  multiples  of  the  central  frequency 
of  the  input;  the  (ideal)  low-pass  filter  eliminates  all  these  zones  except  the 
low-frequency  one,  corresponding  to  the  zeroth  harmonic.  A block  diagram 
is  shown  in  Fig.  3,  together  with  typical  waveforms  and  spectra  at  various 
points  in  the  circuit. 


The  system  treated  here  is  of  course,  an  idealization  of  actual  ones. 
Two  of  the  more  critical  approximations  are  those  of  representing  actual 
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modulations,  which  are  functions  of  time,  by  a constant  value  and  represent- 
ing the  dynamic  detector  characteristic  by  a simple  power  law,  without  in- 
cluding saturation  effects  at  high  voltages  and  smoothing  of  the  discontinuity 
in  the  slope  at  the  origin.  Nevertheless,  the  idealized  system  in  most  appli- 
cations is  a sufficiently  close  representation  to  warrant  analysis. 

Since  there  is  noise  present,  a complete  analysis  should  give  the  prob- 
ability distributions  of  the  output.  This  is  a difficult  problem  which  has  been 
solved  only  for  the  quadratic  detector  { v = 2),  [96,97]  an  important  and  par- 
ticularly simple  special  case.  We  shall  calculate  the  correlation  function  of 
the  output,  a much  simpler  problem  and  one  whose  solution  yields  considerable 
insight  into  the  effects  of  noise  in  the  system.  Two  different  methods  can  be 
applied,  corresponding  to  the  different  statistical  descriptions  of  the  input  by 
the  distributions  of  either  the  instantaneous  value  or  of  the  envelope  and  phase. 
The  former  of  these  is  better  suited  to  the  major  problem,  but  the  envelope 
and  phase  description  give  a better  picture  of  the  zonal  structure,  so  we 
shall  use  it  initially. 


4.2  Output  power 

The  output  can  be  expanded  in  a Fourier  series, 

oo 


I = g[Rcos(<*)t  - 0)]  = / € a (R)cosn(wt  - 0),  (4.2) 

o /-  \ n la  o 


a„(R)  . 


^ /■' 

2tt  / 


n=0 


g [R  cos  (w  t - 0)]  cos  n(u  t - 0)d  u t , 


in  which  c is  the  Neumann  discontinuous  factor,  c = 1,  c =2  for  n > 1. 
n o n — 

The  envelope  and  phase  are  also  functions  of  time,  but  they  vary  so  little 
during  one  period  of  cos  o^t  that  they  may  be  considered  constant  in  calcu- 
lating a^(R). 

The  output  correlation  function  is  then 


TR189 


-37- 
oo 

R(t)  = <I^l2>  av  = ‘n  < av' 

n = 0 

(4.  3) 

Note  that  only  the  phase  difference  occurs  in  this  equation,  so  that  the  re- 
quired probability  distribution  can  always  be  found,  as  in  Eq.  (2.22),  from 
the  distribution  of  instantaneous  values.  The  power  in  each  spectral  zone, 

P , can  be  found  from  the  correlation  function  for  t = 0,P=e  /a  (R)  > 
n n n \ n / av 

since  now  0^  = 0j^.  The  d-c  power  equals  the  square  of  the  mean  value. 

P,  = <Ta  (R)>  . For  the  half-wave  detector,  we  find 

d-c  ^o'  av 

a 
n 

Equation  (4.4)  shows  that  the  detector  output  possesses  the  same  dependence 
on  the  input  envelope  for  all  zones,  so  that  the  only  difference  between  zones , 
is  in  a multiplying  factor,  Blackman  [98]  has  previously  pointed  out  this 
result,  though  stating  the  implication  in  less  general  terms.  This  is  that 
the  relative  power  per  zone  is  a constant  independent  of  the  statistics,  i.e.  , 
regardless  of  what  the  input  is,  the  detector  places  the  same  percentage  of 
the  total  output  power  in  each  spectral  zone. 

Such  a result  is  worthy  of  note,  since  it  is  by  no  means  intuitively 
obvious.  Indeed,  it  seems  possible  that  some  input  could  be  found  that 
would  contribute  an  extra  share  of  power  to  the  higher  zones,  and  a different 
input  whose  power  would  be  distributed  mostly  to  the  lower  zones.  For  the 
V th-law  detector,  no  such  input  existsihowever,  this  detector  is  apparently 
a very  special  case . * 


♦Clearly,  a sufficient  condition  that  the  result  should  hold  is  that  g(Rcos  G)  = 
g^(R)g2(cos  0),  and  the  v-th  law  detector  does,  in  fact,  satisfy  this.  Neces- 
sary and  sufficient  conditions  are  not  obvious;  however,  it  seems  likely  that 
they  are  somewhat  but  not  very  much  weaker  than  the  above,  the  point  being 
that  if  the  envelope  and  cosine  factors  can  not  be  separate  before  the  integ- 
ration is  performed,  n and  R are  likely  to  be  inextricably  mingled  in  the  re- 
sult. As  an  example,  if  a half- wave  rectifier  with  saturation  is  represented  by 
g(v)  = p[  1 -exp(  - av)]  , then  a^  (R)  = p[  I (a  R)]  , and  the  relative  power  in  the  zones 
quite  definitely  depends  on  the  input^slaflstics . 


PR  P(  1) 
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The  low-frequency  zone  is  the  most  important  in  the  treating  of  the 
receiver  problem,  so  that  all  results  are  given  in  terms  of  it.  These  re- 
sults can  be  altered  to  give  the  power  in  a different  zone,  if  desired,  by 
multiplication  by  an  appropriate  constant  factor.  The  averages  are  evalu- 
ated in  Appendix  V;  for  the  total  low-frequency  power,  we  have 


p = 1.)  [ F (-v;l;-p)-f-^°'^.^J  --  F (-v+2;l;- 

° + 1)  ^ ^ ^ 3!  ^ ^ 


4r^(%  + 1) 

and,  for  the  d-c  power, 


p)+o(y"^ 


(4.5) 


‘d-c 


- ^Y‘ 


Iv+l) 


4r^ 


+ 1) 


(ir 


■40^3  - 


^F^(-^;l;-p)  + 


-2. 


^F^(--j+2;l;-p)  +0{  Y ) 


(4.6) 

where  p = A^/2  4iis  the  input  signal-to-noise  ratio.  If  v in  idq.  (4.5)  or 
v/2  in  Eq.  (4.6)  is  a positive  integer,  the  nearly  gaussian  series  termin- 
ates and  the  above  expressions  become  exact;  furthermore,  the  confluent 
hypergeometric  functions  reduce  to  polynomials. 


These  equations  are  plotted  in  Figs.  4 and  5 for  gaussian  noise  and 
for  a nearly  gaussian  noise  of  comparatively  low  density.  As  these  graphs 
show,  the  change  in  the  output  power  caused  by  the  nongaussian  statistics 
is  small  for  small  values  of  v , and  small  for  large  values  of  p.  When  v 
is  greater  than  one,  the  nongaussian  output  power,  (Fig.  4),  is  greater 
than  in  the  gaussian  case;  while  for  v less  than  one,  either  nuay  predomi- 
nate, depending  on  the  signal-to-noise  ratio. 


When  V and  p are  both  less  than  one,  the  nongaussian  power  is  always 
less  than  the  corresponding  gaussian  one;  however,  if  p is  greater  than  one, 
the  nongaussian  power  is  first  greater,  then  less,  the  crossover  point  in- 
creasing from  V equals  zero  to  v equals  one  as  p increases  from  one  to  in- 
finity. For  large  signal-to-noise  ratios,  the  output  power  is  composed 
mostly  of  S X S (signal  X signal)  and  S X N (signal  X noise)  interaction,  while 
with  small  signal-to-noise  ratios,  the  power  arises  mostly  from  n X n 
(noise  X noise)  interaction.  There  is  no  a priori  reason  to  believe  that  the 


FIG.  4 TOTAL  LOW-FREQUENCY  OUTPUT  POWER 
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nongaussian  result  will  bear  the  same  relation  to  the  gaussian  in  the  one 
case  as  in  the  other,  and  in  fact,  Eq.  (4.5)  shows  that  they  do  not.  Ana- 
lytically, these  facts  depend  on  the  behavior  of  the  confluent  byhypergeo- 
metric  functions,  which  are  discussed  in  more  detail  in  Appendix  V.  The 
discussion  on  the  total  power  carries  over  immediately  to  the  d-c  power 
except  that  all  the  values  of  v previously  mentioned  must  now  be  multi- 
plied by  two. 

The  asymptotic  series  for  the  hypergeometric  functions  of  Eqs.  (4.5) 
and  (4.  6)  show  that  the  leading  nongaussian  term  is  of  order  p while 
the  leading  noise  term  is  of  order  p \ Thus,  in  the  strong  signal  case, 
as  well  as  the  usual  suppression  of  noise  relative  to  signal,  there  is  an 
additional  suppression  of  deviations  from  gaussian  statistics  (except  for 
values  of  v near  1).  The  figures  and  equations  show  that  in  order  to  detect 
the  type  of  statistics  present  in  an  experimental  problem  by  a power  meas- 
urement, the  detector  used  should  have  a high  power  law  and  the  noise  power 
should  be  measured  without  any  signal  present.  The  first  of  these  require- 
ments is  intuitively  reasonable,  because,  since  the  tails  of  a distribution 
are  particularly  sensitive  to  the  exact  type  of  statistics,  a large  value  of 

V , which  accentuates  the  effect  of  the  tails,  is  desirable. 

From  a practical  point  of  view,  however,  large  values  of  v will  be 
more  sensitive  to  saturation  effects.  In  an  actual  detector  the  maximum 
output  will  always  be  limited,  either  from  saturation  in  the  detector  itself 
or  in  the  i-f  amplifier  preceding  it.  The  law  of  the  detector,  of  course, 
has  no  bearing  on  the  overloading  of  the  i-f  amplifier,  but  large  values  of 

V will  require  the  detector  to  produce  large  outputs  without  distortion  when 
only  moderate  inputs  are  present.  Four  or  five  would  probably  be  as  high 
a value  of  i'  as  one  would  want  to  choose. 

The  second  requirement  above  is  not  obvious,  since  a locally  injected 
signal  of  known  strength  increase  the  output  noise  power  because  of  the 
cross-modulation  products  of  the  signal  and  the  noise  present  in  the  output 
of  the  detector.  For  low  input  noise  levels,  it  may  be  important  to  be  able 
to  adjust  the  output  nruse  power  to  a level  large  enough  to  be  convenient. 

One  might  think  that  cross-modulation  would  also  accentuate  the  difference 
between  gaussian  and  nongaussian  noise;  however,  the  above  figures  and 
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equations  show  that  quite  the  contrary  is  the  case;  in  fact,  the  power  differ- 
ence is  not  only  less  relative  to  the  total  output  power  with  increasing  signal- 
to-noise  ratio,  but  even  decreases  in  numerical  value.  Thus,  although  adding 
a signal  may  be  worth while  in  order  to  measure  noise  power,  the  method  is 
entirely  unsuitable  for  distinguishing  different  types  of  statistics. 


4.  3 The  output  correlation  function 


For  evaluat  ng  the  correlation  function,  the  distribution  of  the  instan- 
taneous value  is  preferable.  Whila  the  envelope  and  phase  distributions  can 
be  used,  the  final  result  is  then  in  gi  less  tractable  form.  Thr  correlation 
function  can  be  expressed  in  terms  of  the  characteristic  function  by  using  the 
complex  Fourier  transform  of  the  dynamic  characte'ristic  [99,100]  , 


= / g{V)e"^^^dV, 


-oo 


g(v) 


(4.7) 


where  g(V)  is  zero  for  V less  than  some  V^,  and  of  no  greater  than  ex- 
ponential order  at  infinity,  and  C is  a horizontal  contour  in  the  complex  z- 
plane  lying  below  the  singularities  of  f(iz).  The  output  correlation  further 
becomes 


= <g{V^)g(V2)> 


<exp(iziVj  + i2^V2)X^ 


(4.8) 


and  the  average  of  the  exponential  defines  the  second-order  characteristic 
function  (Eq.  2.2). 
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The  input  signal  to  the  detector  in  this  case  is  the  sum  of  a sine  wave  and 
random  noise  which  are  independent  of  each  other,  so  that  the  characteristic 
function  in  Eq.  (4.8)  is  the  product  of  the  separate  characteristic  functions 
of  the  signal  and  the  noise.  Since  both  inputs  are  narrow-band,  the  two  char- 
acteristic functions  are  expanded  in  Fourier  series  in  order  to  exhibit  the 
zonal  structure  of  the  correlation  function;  for  the  detector  problem,  we 
need  only  the  correlation  function  of  the  low-frequency  zone,  R(t). 

The  detailed  evaluation  is  carried  out  in  Appendix  VI;  here  we  men- 
tion that  two  forms  of  the  correlation  function  are  possible,  a convergent 
power  series  in  p,  the  input  signal-to-noise  power  ratio,  useful  when  the 
signal  is  small,  and  an  asymptotic  series  in  1 /p  useful  when  the  signal  is  large 
relative  to  the  noise.  These  series  are  given  in  Eq.  (VI.  13)  and  (VI.  14)  of 
the  appendix,  respectively.  The  general  results  are  sufficiently  complicated 
so  that  a study  of  various  special  cases  is  deqirable. 

The  leading  terms  of  the  series  in  p can  easily  be  found  from  Eq.  (VI. 15); 
however,  the  equation  remains  difficult  to  interpret  unless  the  hypergeometric 
functions  involved  reduce  to  simpler  functions.  The  two  most  important 
values  of  v are  one  and  two,  corresponding  to  the  linear  and  square-law 
detector,  respectively.  The  square-law  case  is  particularly  simple,  viz., 


.2,2 

4^ 


R(t)  = 1 + ^-o  + 


+ 2p(l+  r^)  + P ] ; 


V = 2. 


(4.9) 


This  result  is  exact  for  all  p and  y,  all  other  terms  in  the  general  result 
vanishing.  For  this  particular  value  of  v , the  output  voltage  is  the  same, 
except  for  a scale  factor,  as  that  from  a full-wave  squaring  device;  ac- 
cordingly, the  correlation  function  can  be  calculated  directly  as  a moment. 


2 2 2 

R„(t)  = ^ V(t  ) + A cos  w t ] [ V(t  + t)  + A cos  w (t  + t)]  ^ 

T 4.^  0 o o o-*  ' o o o o •*  '^av 

(4.10) 

This  is  by  far  the  siraplest  way  to  obtain  the  output  correlation  function; 
unfortunately,  the  method  can  only  be  applied  when  v is  an  even  integer. 


For  the  linear  detector,  the  hypergeometric  functions  reduce  to  complete 
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elliptic  integrals,  and  the  leading  terms  are 


R(.)  = ^ 

ir 


40 


E + r B (K-2E)  - 


A,,(.) 


-ai 2r„D+  , (2D  +K) 

124/  ° 12+'^ 


+ P 


“ 12?^ 


E ^ 2E-B 
4E  - 2K 2“  + "’o  •; TT 


1-r 


1-r 


124<^ 


Er  - B 
o 

T7“ 


A22<'> 


I24J 


E + 2B 

-j-  + Tq[ y-  2C] 


1-. 


1-r 


+ 0{y'^.p^)  \ 


v=l. 


(4.11) 


In  this  equation,  the  functional  dependence  of  r^{t)  on  t,  and  of  the 
elliptic  integrals  on  r^,  e.g.,  E(r^),  has  not  been  explicitly  expressed. 

If  it  is  thought  preferable,  the  various  elliptic  integrals  appearing  can 
be  expressed  in  terms  of  any  two,  preferably  K (r^)  and  E (r^)  the  complete 
elliptic  integrals  of  the  first  and  second  kind  respectively,  by  using  the  re- 
lations [ lOl] 


E + (l-r^)K  = r^B, 
o o 

K - E = r^D, 
o 

(2-r^)  K = 2E  + r^C.  (4.12) 

0 0 

All  five  integrals  have  been  tabulated,  e.g.,  in  Jahnke  and  Emde . This 
result  for  the  linear  detector  is  not  exact,  either  in  p or  in  y;  furthermore, 
there  is  no  essentially  simpler  way  of  obtaining  this  result,  in  contrast  to 
the  case  of  the  quadratic  detector. 


The  correlation  function  when  the  signal  is  much  stronger  than  the 
noise  is  considerably  simpler  than  for  the  v/eak-signal  case.  The  leading 
terms  of  the  result  for  a vth-law  detector  are 


R{t)^ 


p2r2{ 


v+1) 


4r^^+i) 


P^ 


T 1 + r 


% 


t 
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JL 


^ V , V 


i)(A^i)  + (A^^)  + 2(j)^  + {j  - i)^J 


+ 0(p"^,  y-h 


(4.1?) 


The  correlation  function  shows  the  same  suppression  effects  previously- 
mentioned  in  connection  with  the  output  power,  viz.  , the  first  noise  term 
i&  of  lower  order  (in  p)  than  the  signal,  and  the  first  non-gaussian  term 
is  of  lower  order  than  the  first  gaussian  one. 


4 .4  Noise  models 

To  proceed  further  requires  explicit  models  of  the  noise  statistics 
in  order  to  fix  the  time  dependence  of  the  semi-invariants.  An  immense 
variety  of  types  is  naturally  possible;  three  models  among  those  felt  to  be 
the  most  important  are  treated  here.  The  amplitude,  shape  of  the  envelope, 
and  the  phase  of  the  individual  pulses  are  able  to  affect  the  time  dependence 
of  the  statistics,  while,  in  our  models,  time  dependence  enters  through  the 
shape  of  the  pulse  envelope  alone.  This  apparent  restriction  is  justified  as 
being  usually  true  of  actual  noise  sources  (cf.  Appendix  IV  where  this  point 
is  treated  in  connection  with  radar  clutter)  and  as  being  irrelevant  to  the 
real  problem  of  choosing  a suitable  time  variation.  By  this  last  is  meant 
that  the  time  dependence  of  the  semi-invariants  is  what  really  matters  and 
this  can  remain  the  same  whether  ascribed  to  pulse  envelope,  amplitude, 
phase,  or  some  combination  of  all  three. 

Pertinent  data  on  the  three  models  are  summarized  in  Table  4.  1.  The 
exponential  pulses  represent  impulses  after  having  passed  through  a single- 
tuned  circuit,  as  might  occur  with  precipitation  static  in  an  aircraft  receiver 
of  shrimp  noise  in  a sonar. 

Although  a single  tuned  circuit  is  not  an  accurate  model  of  the  tuned 
stages  of  a receiver,  this  type  cf  time  dependence  is  important  because  it  is 
necessary  (and  sufficient)  if  the  process  is  to  be  Markoffian  in  the  limit  of 
increasing  density  [ 102]  . 

The  pass  band  of  an  actual  receiver  is  an  involved  function  of  the 
number  of  stages  and  the  exact  coupling  network  between  stages.  As 
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an  abstraction  from  the  details  associated  with  any  particular  i-f  strip, 
we  shall  take  a pass  band  of  gaussian  shape  as  our  model,  which  preserves 
the  essential  features  of  a pass  band  while  remaining  analytically  tractable. 
Admittedly,  a gaussian  pass  band  is  not  physically  realizable,  but  it  is  a 
good  approximation  to  the  magnitude -frequency  curve  of  actual  amplifiers 
(if  not  to  the  phase-frequency  curve),  and  possesses  the  important  virtue  of 
simplicity.*  Impulses  passing  through  this  i-f  amplifier  will  become 
gaussian  pulses  as  in  the  second  model  of  Table  4.  1. 

The  third  type  chosen  is  that  of  square:  pulse  envelopes  of  finite  duration. 
This  illustrates  noise  whose  values  are  independent  when  separated  by  a 
sufficiently  long,  but  finite  time.  The  model  fits  the  sonar  and  radar  clutter 
problem  when  relative  motion  between  the  transceiver  and  scatterers  is  slight. 

Notice  that,  for  these  three  cases,  all  the  semi-invariants  can  be  ex- 
pressed as  functions,  in  fact  powers,  of  the  input  correlation.  Accordingly, 
the  time  enters  only  on  a normalized  scale  and  only  implicitly  through  the 
input  correlation  function.  Figure  6 shows  the  time  dependence  of  the  fourth- 
order  semi-invariants  compared  to  the  correlation  function,  showing  the  ex- 
tent of  the  more  rapid  decrease  to  zero  of  the  higher  semi-invariants.**  The 
corresponding  figure  for  the  linear  model  is  obvious. 


4.5  Results  for  the  quadratic  and  linear  detectors 

Table  4.Z  lists  the  correlation  functions  after  a quadratic  detector  for 
the  three  different  noise  models.  The  unnormalized  values  of  the  functions 
in  the  table  are  not  of  primary  interest  here  since  the  parallel  problem  con- 
cerning output  power  has  been  discussed  earlier.  The  correlation  function 


*H.  Wallman  has  shown  that  the  pass  band  of  cascaded  networks  whose  individual 
step  function  responses  have  no  overshoot  tends  toward  the  gaus  sian  (Second  Sym- 
posiumin  Applied  Mathematics,  Am.  Math.  Soc . 1950,  p.  91).  Furthermore  it 
seems  likely  that  one  can  extend  these  results  to  any  cascaded  network  except 
to  those  tuned  for  a Butterworth  response. 

**The  linear  model  is  one  extreme  of  possible  time  behavior  in  that  its  semi- 
invariants of  all  orders  decrease  no  faster  than  the  correlation  function.  One 
may  conjecture  that  no -A.22(t)  can  decrease  faster  than  the  square  of  the  correla- 
tion function,  in  which  case  the  other  two  models  represent  the  other  extreme ; 
however,  we  have  been  unable  to  prove  this. 


NOISE  MODELS 


2 ,m+n 


TABLE  4.2 


CORRELATION  FUNCTIONS  FOR 
THE  QUADRATIC  DETECTOR 


1.  Exponential  Model:  = e 


-P|t| 


.2,2 


R(t)  = ( ( 1 + p)2  + 2pr^^  + ( 1 + i)r^] 


2.  Gaussian  Model: 

,2.2 


r s e 
o 


■~T~ 


R(t)  = [ ( 1 + p)2  + 2p  + ( 1 + — )r  J] 


3.  Linear  Model:  = |l-  Plt|  , p |t|  1 

1 0 .Pltl^l 

2 2 

R(t)  = [ ( 1 + P)^  + (2p  + ^)  + r J ] 
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also  gives  information  about  statistical  dependence  in  time  (or  equivalently, 
the  spectral  distribution  of  power)  which  is  most  clearly  exhibited  by  normal- 
izing thecorrelation  function  so  that  it  varies  from  one  to  zero.  Accordingly, 
we  define 


r 


= «-(t)  - R(qo) 
out'  ' R(0)  - R(oo) 


(4.14) 


The  more  rapidly  r^^^(t)  drops  to  zero,  the  less  correlation  there  is 
between  successive  samples  of  the  output  wave. 


Some  statements  about  the  general  effects  of  signal-to-noise  ratio  and 
noise  density  can  be  made  immediately.  Passing  through  a nonlinear  device 
cannot  increase  the  correlation  present  in  the  input  wave  , whence  ^Qut^^^ 
regardless  of  the  values  of  p or  y.  Furthermore,  the  noise  suppression  effect 
when  p is  large  means  that  the  most  rapid  decrease  of  'wiU  occur  when 

no  signal  is  present  and  that,  with  increasing  signal  ^Q^^^(t)  will  tend  to  r^(t)  • 
The  nongaussian  suppression  effect  means  that  the  difference  between 
for  gaussian  and  nongaussian  noise  of  fixed  density  and  equal  input  power  will 
decrease  as  the  signal  is  increased.  The  effect  of  nongaussian  noise  on  the 
amount  of  correlation  in  the  output  depends  in  an  involved  fashion  on  the 
detector  law,  signal-to-noise  ratio,  and  noise  model,  so  that  the  discussion  of 
this  point  is  best  deferred  till  after  the  presentation  of  some  specific  results. 


Figures  7,8  and  9 show  r^^^(t)  for  a quadratic  detector  and  the  three 

different  noise  models.  Note  that  the  three  general  features  of  r (t)  men- 

tioned  above  are  borne  out  by  the  figures.  For  the  quadratic  detector,  the 

exponential  and  gaussian  models  give  less  correlation  with  nongaussian  noise 

than  with  gaussian,  while  for  the  linear  model,  the  opposite  is  true.  Because 

of  the  simplicity  of  the  expression  for  the  correlation  function,  we  can  see 

quite  easily  why  this  is  so.  Gaussian  noise  after  the  quadratic  detector  has 

one  undistorted  noise  term  (2pr  , arising  from  s x n intermodulation)  and 

2 ° 

one  scrambled  noise  term  from  n x n modulation).  If  the  noise  is  non- 

gaussian an  additional  n x n term  enters,  containing  the  semi -invariant^22(^)  • 
For  the  exponential  and  gaussian  models,  this  term  adds  to  the  scrambled 
noise  to  make  the  nongaussian  result  less  correlated  than  the  corresponding 
gaussian  one,  while  for  the  linear  model,  the  semi -invariant  term  adds  to  the 
undistorted  noise  term,  acting  to  increase  the  correlation  over  the  gaussian 
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arnount. 

Table  4.3  shows  the  correlation  functions  for  the  linear  detector.  Be- 
cause the  expressions  do  not  terminate,  it  is  not  possible  to  cover  the  entire 
range  of  values  of  p and  *y  . Strongly  nongaussian  noise  (7=  1)  can  be  included., 
however,  because  of  the  numerical  factors  involved  in  the  Edgeworth  series. 

The  first  correction  term  is  about  the  size  of  2 7V.^q/4!4»  or  l/4y.  The  second 
is  in  two  parts,  one  half  the  square  of  the  first  term  or  and  one  involving 

33  2 

sixth-order  semi -invariants , 2^A^q/6!4<  or  l/36\  , both  of  which  are  suf- 
ficiently small  to  be  neglected  inasmuch  as  the  first  correction  is  itself  not 
large  compared  to  the  gaussian  part  of  the  results.  The  correlation  function 
is  much  more  sensitive  to  the  dependence  on  signal -to -noise  ratio.  The  strong 
signal  series  is  sufficiently  accurate  when  p is  greater  than  two;  the  weak  signal 
series,  however,  can  only  be  used  for  p = 0.5  at  the  most. 

The  normalized  output  correlation  functions  are  shown  for  the  linear 
detector  in  Figs.  10-12.  Previously,  with  the  quadratic  detector  the  nongaussian 
correlation  function  was  either  always  greater  or  always  less  than  the  correspond- 
ing gaussian  one.  For  the  linear  detector,  however,  the  correlation  is  greater 
for  nongaussian  noise  when  the  signal  is  weak  and  greater  for  gaussian  noise 
when  the  signal  is  strong.  An  additional  complication  is  that  the  curves  for 
gaussian  and  nongaussian  noise  cross  each  other  near  the  origin  of  time  when 
p is  large,  although  for  by  far  the  most  of  the  time  the  nongaussian  curve  lies 
below.  The  curves  again  illustrate  that  the  difference  between  the  two  types 
of  noise  is  greatest  when  the  signal  is  small,  and  that  the  difference  is  sup- 
pressed when  a strong  signal  is  present.  For  the  linear  model  the  suppression 
is  of  higher  order,  so  that  the  difference  is  particularly  small. 

4 . 6 The  effect  of  nongaussian  statistics  in  general 

The  effect  of  the  nongaussian  statistics  on  varied  in  these 

two  cases,  depending  in  a rather  obscure  manner  on  the  law  of  the  detector, 
noise  model  and  signal -to -noise  ratio,  so  that  an  investigation  of  general 
values  of  v is  worth  while;  however,  since  the  linear  and  quadratic  detectors 
are  by  far  the  most  important,  a qualitative  look  at  the  results  for  arbitrary 
values  of  v suffices  to  complete  the  discussion  of  the  output  cor  relation  function. 
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The  primary  point  in  question  is  whether  nongaussian  noise  increases  or 
decreases  the  correlation  of  the  output  compared  to  gaussian  noise  of  the 
same  input  power.  We  have  obtained  the  sign  of  the  change  in  correlation 
as  a function  of  v and  p for  the  special  cases  above  and  for  various  limiting 
values  of  v and  p,  but  have  not  been  able  to  get  a really  satisfactory  assurance 
that  no  additional  changes  of  sign  occur  for  intermediate  values  of  v or  p. 

The  best  way  to  present  the  result  is  to  show  the  first  quadrant  of  a(v,p)- 
plane  in  which  the  boundaries  and  regions  of  increased  or  decreased  cor- 
relation are  explicitly  labelled. 

The  behavior  along  the  lines  v = 1 and  v = 2 is  known  from  the  results 
of  Figs.  7-12;  and,  in  addition,  we  have  found  the  sign  of  the  change  in  cor- 
relation for  large  p and  v,  and  for  p = 0 and  all  v.  Appendix  VII  gives  the 

explicit  forms  of  r ^(t)  from  which  Fig.  13  was  constructed. 

^ out  ® 

The  dependence  of  the  change  in  the  total  output  power  and  in  the  d-c 
power  on  the  law  of  the  detector  leads  one  to  expect  that  the  normalized  out- 
put correlation  function  should  also  exhibit  some  such  settled  type  of  be- 
havior when  V is  larger  than  some  number,  perhaps  not  exactly  two.  In  fact, 
Fig.  11  bears  this  out.  Nonetheless,  while  all  the  information  that  appears 
in  the  figure  is  correct,  and  any  additional  changes  in  the  behavior  of  r^^^(t) 
seem  unlikely,  it  cannot  be  categorically  stated  that  all  the  changes  in  cor- 
relation are  included. 

The  figure  shows  that,  with  the  exponential  and  gaussian  models,  when 
the  law  of  the  detector  is  greater  than  two,  the  correlation  of  the  output  noise 
is  greater  when  the  noise  is  nongaussian.  Only  when  the  detector  law  is  less 
than  two  does  the  relative  strength  of  the  signal  affect  the  result.  Near  v = 1, 
in  the  shaded  area,  the  output  correlation  functions  for  gaussian  and  non- 
gaussian noise  cross  each  other  so  no  definite  assignment  of  this  region  can 
be  made.  Figures  10-12  for  the  linear  detector  show,  however,  that  for  most 
of  the  time,  the  correlation  of  nongaussian  noise  is  less  than  that  of  gaussian. 
When  V is  less  than  two,  the  nongaussian  noise  terms  have  signs  which  differ, 
giving  a noise  cancellation  effect  with  strong  and  weak  signals  (for  the  specific 
equations,  see  Appendix  VII). 
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The  effect  is  not  exactly  the  same  in  the  two  cases  becavise  the  noise  is 
different,  coming  mostly  from  signal-by-noise  modulation  products  when 
the  signal  is  strong  and  from  noise-by -noise  modulation  products  when  the 
noi se  is  weak. 

The  linear  model  differs  from  the  other  two  in  the  decrease  of  the  semi- 
invariants with  time  compared  to  the  input  correlation  function.  Just  as  in  the 
previous  discussion  of  the  linear  and  quadratic  detectors,  one  expects  the  cor- 
relation of  this  model  to  be  greater  than  that  of  the  other  two.  Figure  13  shows 
that  here  the  nongaussian  noise  is  more  strongly  correlated  for  all  values  of 
V when  no  signal  is  present,  and  that  the  band  of  increased  correlation  at  high 
signal  levels  has  moved  to  larger  values  of  v . 

The  magnitude  of  the  change  in  correlation  will  be  greater  when  the 
signal  is  weak  or  absent.  At  large  signal  levels,  the  change  is  suppressed 
by  a factor  of  1/p  for  the  exponential  and  gaussian  models  and  by  a factor 
of  1/p^  for  the  linear  model. 

4.  7 The  smoothing  effect  of  finite  time  averages 

Up  to  now,  the  discussion  of  the  output  of  the  detector  has  been  con- 
cerned with  average  values,  which  implies  averaging  over  all  time.  When 
the  output  is  averaged  over  a finite  time,  the  result  is  not  a parameter  of 
the  output  distribution  as  for  the  infinite  average,  but  rather  again  a random 
variable,  although  one  with  a more  concentrated  distribution  than  before. 

In  actual  practice,  of  course,  only  a finite  time  is  available  in  which  to  operate 
on  the  output  of  the  detector.  The  finite  averaging  time  makes  an  analysis 
more  difficult,  but  the  variance  of  the  averaged  output  can  be  found  [ 103]  . 

The  effect  of  a finite  averaging  time  depends  partly  on  the  video  filter. 
This  is  a low-pass  filter  transmitting  the  low-frequency  spectral  zone  and 
cutting  off  the  spectral  zones  at  the  carrier  frequency  and  its  multiples,  and 
is  followed  in  turn  by  an  ideal  finite  time  integrator.  From  the  reference 
quoted  above,  we  have  the  result  that  the  variable  z = (1/T)  J*q  lQ(t)dt, 
where  T is  the  averaging  time  and  is  the  low-frequency  output  of  the  detector, 
has  as  its  variance. 


EXPONENTIAL  AND  GAUSSIAN 
NOISE  MODELS 


0 2 1/  4 6 


fig.  13  EFFECT  OFNON  GAUSSIAN  STATISTICS  ON  THE 
OUTPUT  CORRELATION 
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X 

= [R(t)  -1^]  [1  - tn]  dt 


(4.15) 


where  the  averages  under  the  integral  sign  are  ensemble  averages.  Notice 
that  to  find  the  variance  of  Z requires  knowing  a more  complicated  statistic 
of  the  infinite  averaging  problem. 

The  general  effect  of  the  averaging  can  be  readily  seen.  When  T is 

shorter  than  the  time  in  which  the  output  correlation  function  goes  to  its 

final  value,  the  averaging  is  almost  entirely  ineffective  so  that  the  variance 

of  Z is  as  great  as  that  of  I^.  On  the  other  hand  when  T is  much  greater  than 

the  time  in  which  the  correlation  function  varies,  the  mean  and  variance  of 

Z are  effectively  composed  of  the  sum  of  independent  random  variables , whose 

number  is  proportional  to  the  averaging  time.  Thus  for  long  averagingtim.es, 

2 2 

the  relative  variance  cr  fl  , will  vanish  as  l/T . Explicit  results  are  in  general 
difficult  to  obtain  because  t appears  implicitly  in  the  correlation  function; 
however,  for  the  quadratic  detector,  the  elementary  form  of  the  output  cor- 
relation function  again  leads  to  particularly  simple  results. 

The  relative  variance,  cr^ / I is  the  best  measure  of  the  elimination 
of  noise  effected  by  averaging.  Table  4.4  gives  the  quadratic  detector  equa- 
tions, and  Figs.  14,  15,  and  16  present  graphs  of  the  variation  of 
the  variance  with  time.  Just  as  before,  with  increasing  signal,  the  relative 
noise  output  becomes  less  and  the  difference  between  gaussian  and nongaussian 
noise  also  is  less  pronounced.  However,  this  orderly  sequence  is  disturbed 
for  weak  signals  and  long  smoothing  times.  The  smoothing  is  more  effective 
against  the  n x n component  of  the  relative  variance  than  the  s x n component 
because  the  n x n part  is  less  correlated  than  the  other;  thus,  with  long  aver- 
aging times,  the  remaining  fluctuation  is  mostly  from  the  a x n component. 

At  low  signal  levels,  where  the  variance,  even  with  no  smoothing,  is  only 
slightly  less  than  that  with  no  signal,  the  variance  with  long  smoothing  times 
will  be  greater  because  the  "no-signal”  noise,  arising  entirely  from  n x n 
interaction,  is  more  strongly  suppressed. 

Another  interesting  point  is  that  the  nearly  gaussian  noise  always  has 
a larger  variance  tnan  gaussian  for  all  three  models,  when  noise  alone  is 
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present.  This  behavior  is  in  contrast  to  that  of  the  correlation  function 
where  the  linear  model  differs  from  the  other  two.  In  fact,  for  finite  aver- 
aging,  the  three  models  behave  the  same  nearly  everywhere  in  the  (v,p)- 
plane.  Appendix  VIII  takes  up  the  question  in  more  detail;  the  results  are 
that  nearly  gaussian  noise  has  a greater  variance  than  gaussian  noise  of  the 
same  input  mean  square,  except  when  the  signal  is  large  and  v lies  in  the 
approximate  range  0.5  to  1.5,  in  which  case  the  nearly  gaussian  variance 
is  less.  The  boundary  of  the  two  regions  is  not  quite  the  same  for  the 
three  different  models,  but  this  is  the  only  difference  between  them. 


[(2p+M  (3pr  - 1)  + j{4^T  - 1)]  ,pT 


FIGURE  14 


=1.565 


FIGURE  15 


FIGURE  16 
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V 

Conclusion 

In  the  preceding  sections,  we  have  discussed  the  form  of  the  nearly 
gaussian  distributions  and  the  physical  noise  sources  that  may  be  expected 
to  produce  nearly  normal  random  noise.  In  the  previous  section,  the  problem 
of  the  rectification  of  nongaussian  noise  has  been  considered  and  most  of  the 
important  statistical  parameters  of  the  output  determined  in  general,  although 
in  perhaps  intractable  form.  The  qualitative  effect  of  nongaussian  statistics 
on  the  output  power,  output  correlation  function,  and  finitely  averaged  output 
have  been  determined  in  reasonable  generality,  with  explicit  formulae  and 
graphs  for  the  most  important  special  cases. 
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Appendix  I 
Stationarity 

Even  though  the  instantaneous  value  of  the  noise  wave  is  stationary, 
in  general  the  slowly  varying  components  are  not;  however,  since  the  noise 
is  narrow-band,  to  a satisfactory  degree  of  approximation,  they  may  be  as- 
sumed so.  The  difficulty  is  the  same  as  with  a modulated  carrier,  viz.  , 
a phase  shift  is  involved  which  spoils  stationarity  unless  the  phase  distribu- 
tion is  uniform. 

The  importance  of  establishing  stationarity  lies  in  the  fact  that  analytic 
results  are  obtained  from  statistical  averages,  while  experimental  results 
are  obtained  from  time  averages.  If  the  random  process  involved  is  stationary 
or,  more  precisely,  is  composed  of  stationary  subsets  that  possess  no  non- 
trivial stationary  subsets,  then  the  two  kinds  of  averages  are  equal  (to  with- 
in a subset  of  measure  zero). 

To  discuss  stationarity,  it  is  most  convenient  to  use  the  ensemble 
representation.*  For  the  Poisson  ensemble,  we  may  use 

V(t;t',v|j)  = R(t-t’)  cos  (w^[t-t']  - 4<)  > (I*  1) 

where  t'  is  uniformly  distributed. 

Equation(I.l)  is  not  a random  function;  it  is  a particular  member  of  the 
ensemble  fixed  by  the  ensemble  parameters  t'  and  ijj  and  as  such  is  a 
determinate  function.  At  a later  time  f,  the  time  of  Eq . (I.  1),  is  re- 
placed by  t + T-  This  generates  a new  ensemble  which  can  be  compared 
to  the  old  one  only  when  the  time  difference  fis  absorbed  into  the  en- 
semble parameters  by  some  appropriate  transformation.  After  the  two 
ensem.bles  have  the  same  time  value,  they  may  be  compared  to  see  if 
they  are  the  same,  which  is  the  stationarity  requirement.  In  this  case, 

T can  be  added  to  t' , 

V(t  +r;  t',  4i)  = V(t;  t'  - -r.  i|>)  . (1.2) 

Since  t'  is  uniformly  distributed  the  translation  by  T will  leave  the  overall 
ensemble  invariant. 

*An  ensemble  is  a set  of  functions  together  with  a probability  distribution 
definea  over  the  set. 
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For  the  in-phase  and  quadrature  components,  a further  restriction 
is  needed.  We  have  specifically 


X(t;t',4))  = R(t-t')  cos  (00^  t'  + i|j) 
Y(t;t',4i)  = R(t-t')  sin  + i|j) 


(1.3) 


and 

X(t  + r;t',4j)  = X(t;t'  -r,  • 


(1.4) 


This  last  equation  indicates  that,  for  the  slowly  varying  components 
to  be  stationary,  the  phase  must  be  purely  random  as  well  as  the  time  t' ; 
however,  the  narrow-band  character  permits  a relaxation  of  this  resultfor 
nearly  gaussian  distributions. 

Since  we  have  a Poisson  ensemble,  the  semi -invariants  of  X and  Y 
are  simply  related  to  the  moments  of  the  individual  pulses. 

J^mn^  <x”^y">=  y <h”"'‘’'^(p[t-t'])cos"^[a)^t'  + i|j]xsin*^[w^t’  + i|j]  > (1.5) 


The  trigonometric  part  canbe  written  as  a Fourier  series  and  each  term 
interpreted  as  the  sine  or  cosine  transform  of  h”^^^((3[  t-t'] ),  evaluated  as 
a multiple  of  the  central  .frequency,  (the  phase  may  be  absorbed  in  t'  since 
t'  is  uniformly  distributed).  The  pulse  envelope  varies  much  more  slowly 
than  the  central  frequency  part,  so  that  its  transform  will  have  hardly  any 
spectral  components  at  the  central  frequency.  This  will  also  be  true  of  the 
transforms  of  the  smaller  powers  of  the  pulse  envelope,  so  that  all  terms 
in  the  Fourier  series  may  be  neglected  except  the  leading  one,  correspond- 
ing to  the  zeroth  harmonic,  which  gives  the  isotropic  result.  Thus,  for 
nearly  gaussian  noise  where  the  higher-order  moments  are  negligible,  the 
distribution  is  effectively  isotropic.  The  conclusion  can  not  be  extended 
to  a narrow-band  distribution  in  general,  because  of  the  different  behavior 
of  the  high-order  moments  . Since  the  energy  in  an  individual  pulse  is 
finite,  the  pulse  envelope  must  either  last  for  a finite  amount  of  time  or 
vanish  sufficiently  strongly  at  infinity.  Taking  powers  of  envelopes  of  the 
decreasing  class  will  increase  the  rate  of  fall  and  broaden  the  spectrum, 
until,  for  large  values  of  m+n,  the  spectral  components  at  harmonics  of  the 
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central  frequency  are  no  longer  negligible  and  the  distribution  thus  not 
isotropic . 

Table  I.  1 illustrates  the  K.bove  remarks,  for  three  different  types  of 
pulse  envelopes:  (i)  exponential  pulses,  characteristic  of  a single  tuned  cir 
cuit,  (ii)  gaussian  pulses,  characteristic  of  a multi-stage  amplifier,  and 
(iii)  square  pulses  of  finite  duration.  The  duration  of  the  major  part  of  the 
pulse  is  1/p,  so  that  p is  roughly  half  the  bandwidth  of  the  i-f  circuits  in- 
volved; accordingly,  the  ratio  of  harmonic  term  to  isotropic  term  is  ex- 
pressed in  terms  of  Q,  the  reciprocal  of  the  fractional  bandwidth.  Frac- 
tional bandwidths  in  radar  are  about  0.  1,  so  that  the  moment  for  which  a 
harmonic  term  first  reaches  5%  of  the  isotropic  one  is  the  fifth  moment  for 
exponential  pulses  and  the  thirty -fifth  for  gaussian  pulses.  Square  pulses 
do  not  present  the  same  problem,  since  all  the  moments  behave  the  same 
as  far  as  the  relative  contributions  of  harmonics  are  concerned;  for  the 
assumed  Q,  the  first  harmonic  is  5%  of  the  isotropic  term  and  higher 
harmonics  all  contribute  less. 

The  above  remarks  show  that,  for  nearly  gaussian  Poisson  ensembles 
isotropy  follows  from  the  narrow-band  structure,  so  that  the  in-phase  and 
quadrature  components  are  stationary  under  the  same  restrictions  as 
placed  on  the  instantaneous  value;  this  is  sufficient  for  our  purposes. 
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Table  I.  1 

Envelope  Moments  and  Spectra 
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Appendix  II 

Nearly  Gaussian  Distributions  in  Terms  of  Parametric 

Derivatives 


This  series  is  useful  when  the  noise  ensemble  has  a phasje  uniformly 
distributed,  independent  of  the  envelope.  As  has  been  shown,  nearly  gaussian 
narrow -band  noise  waves  have  a structure  equivalent  to  this.  Thus  the 
random  variable  to  be  considered  is 

Z{t)  = R(t)  cos  [o3^t  - 0{t)]  , {U.  1) 


wher  e 


W{R,0)  fw(R)/2Tr  0<  8 < 2tt  , 

0 elsewhere. 

The  first-order  moments  are  easily  calculated. 
ilT 


, cos’^(o)  t - 0)d8/2iT  , 
av  / ' o 


-p-*  i n+ 1 1 


VttP  1 ) 


(U.2) 


^2u+l  “ ° 


and  oo 

k=0 


(II.  3) 


The  odd  semi -invariants  are  also  zero,  as  may  be  seen  by  comparing 
their  power  series,  Eq . (2.6),  with  the  one  for  the  moments.  We  have 
accordingly 
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1 T-  i'c\  _ \ * (i5)  \ _ 1 r , \ ^ ^Zm  >..Zm-, 

In  F^(^)  / ^ -yt-  ^ 1 


(II.  4) 


Since  only  even  powers  appear  in  the  series  for  the  moments,  the  same  must 
be  true  for  the  semi -invariant  series,  which  implies  thatX^j^^^  = 0-  Conse- 
quently,  r _ ^ 


(II.  5) 


The  exponential  can  be  expanded  in  a power  series  except  for  its 
first  term;  the  leading  terms  will  be 

r X ^ -I 

‘ +^u5)'‘  + ^<iB)‘'  + — ri‘5)®  + o(v'^  = ^ 

^ L 2!(4!) 


(II.  6) 


Clearly  this  equation  can  be  written  in  terms  of  derivatives  with  respect  to 
ijj,  viz . : 

r 7 ^ \ ^ d ^ "*"5 

^Z<3)=  l + + — ^ + . (II. 7) 

^ L H 2!(3!)  I 


The  corresponding  series  for  the  second-order  distribution  is  more  in- 
volved, since  the  narrow-band  structure  appears  explicitly  in  the  moments. 
The  second-order  moments  are 

u.  (t)  = <(cos"^(u)  t -0,)cos’^(w  [t  +t]  - 9,)B  , (II. 8) 

^mn  \ 1 2 /av  \ ' o o 1 ' o'-  o ■*  2 /av  ' ' 

which  gives,  after  expanding  the  second  cosine, 

n 

(t)  = <R"?r">  / <cos""^w  t-e)8in^o)  t-0)>. 

^mn  . 1 2^av^/ \ ' o o -'0  av 


m-f  n-k, 


(u)  t -0,)sin^(u)  t -0,)^o  ,(II.9) 

oo  1 oo  l'0jav'  ' 
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wher  e 

6 = 6^  - . (II.  9) 

The  average  over  0^  will  be  different  from  zero  only  when  both  k and 
m +n-  k are  even;  thus  all  odd-order  moments  again  vanish.  As  before, 
the  same  is  true  of  the  semi -invariants . The  average  over  the  phase 
change,  6,  must  be  postponed  until  the  end  in  order  to  use  the  series  of 
derivatives.  The  first  set  of  terms  in  the  semi-invariant  series,  Eq.  (2.5), 
is 

- + 24,r^5i32^.os(a,^t-e)  + ^2  • (II-  10) 

where  the  subscripts  to  the  i)>'s  have  been  inserted  to  permit  differentia- 
tion of  each  of  the  terms  separately.  In  order  to  obtain  a series  in 
parametric  derivatives,  it  must  be  shown  that  the  appropriate  powers 
of  cos(u^t  - 6)  in  a given  term  of  the  expansion  of  the  exponent 

can  be  found  by  taking  appropriate  derivatives  of  the  expression  in  Eq . 

(II.  10).  Since  k in  Eq.  (II. 9)  is  even,  the  sine  factor  can  be  expressed 
in  terms  of  cosines  and  expanded,  so  that  the  m,  nth  moment  term  con- 
tains 5^3?  cos*^  •)((*)  t-0),  where  m+n  and  j are  even,  and  j<n.  Then  by 
taking  n-j  derivatives  with  respect  to  j/2  derivatives  with  respect 

to  3.nd  (m-n+j)/2  derivatives  with  respect  to3^>  the  explicit  dependence 
on  both  the  variables  and  the  narrow-band  structure  is  suppressed.  This 
is  a total  of  (m+n)/2  derivatives,  just  half  the  number  an  Edgeworth  series 
would  require.  Once  assured  that  the  series  can  be  written  down,  it  is 
easier  to  find  the  terms  directly,  rather  than  to  search  for  a general 
formula . 


For  Poisson  noise,  the  leading  terms  are 


^2(^1 ’^2’^)  ^ 


^40  5l  _ A3i(t)cos(u)^t-0)  3 


1 + -- 41-‘  - ^ 


3! 


Azz(t) 


~1T  " ^ co3(uj^t-0)^^;^^4— ^2 


3 A40  >-  4 


X exp|-  j[  2i|>r^5j52^°®  +2^2^^ 
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0 av  ’ 


(II.  11) 
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giving  a.t  once 


1 + 


_19  r d 

31 


[■ 


d‘ 


8^^  8i|i 


X exp 


J-i 


i[  4ij5/+  24,r^Bi52C08  (w^t  - 0)  + 4.252^1 


>Q  • (11.12) 

■^0  av  ' 

‘l>l  = il>2=^ 


After  the  derivatives  have  been  found,  the  subscripts  are  no  longer 
desirable  since  the  noise  is  stationary.  Eq.  (11.  12)  is  the  natural  generaliza- 
tion of  Eq.  (11.7)  as  can  be  seen  by  expressing  the  differential  operator  in 
matrix  form 


gauss 


>« 

4<1-i1>2-‘^ 


^40^31^22 
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9iii 
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^31^22^13 
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^22^13^40 
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(U.13) 
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Appendix  III 

Envelope  and  Phase  Distributions 

Often  the  distributions  of  the  slowly  varying  components  are  needed. 
It  is  obviously  desirable  to  relate  these  to  the  distributions  of  the  instan- 
taneous value  which  are  usually  previously  known.  Again  it  is  most  con- 
venient to  work  with  the  characteristic  function  and  the  series  for  the 
moments.  The  in-phase  and  quadrature  components  are  defined  from 
Eq.  (2.  19),  V(t)  = X(t)  cos  cj^t  + Y(t)  sin  co^t.  The  first-order  moments 
are  then 

p = cos"'  0'  sin^  0'>  (UI.  1) 

'^mn  ^ / av  ' ' av 

The  characteristic  function  when  written  in  polar  coordinates  is 

oo 

_ , j • \ ' ^mn  .m+n  m . . n , 

„{pcos  0,  psind)  = > — ; — i(ip)  cos  0 sin  0 (III.  2) 

M=0 

m+n»M 

Substituting  Eq.  (III.  1)  in  Eq.  (III.  2),  and  using  the  binomial 
theorem  on  the  terms  of  fixed  value  of  M,  we  find  that 


P’  Psin0)  = 


M=o 


(ip) 


M 


M ! 


<R^cos^(0-  -0)>^^ 


(m . 3) 


If  0 is  independent  of  R and  uniformly  distributed,  the  average  may 
be  carried  out  to  give  the  Mth  moment  of  V;  otherwise,  the  shift  in  the 
angle  by  an  amount  0 may  change  the  value  of  the  average  so  that  it  is 
not  th^.  Mth  moment.  In  the.  former  case,  we  have  the  desired  result 

Fx  Y^P  P (in -4) 


To  obtain  the  distribution  of  the  envelope  alone,  the  restriction 
may  be  relaxed.  In  general,  we  have  that 


W(R,0) 


pd  pd0 


^i  pRcos(0-0) 


F„  y(  pcos  0,  psin  0)  . 
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The  integration  over  0 can  be  performed,  after  which  d will  appear  only 
in  the  characteristic  function.  Then  the  cosine  factor  in  Eq.  (III.  3)  can 
be  averaged  over  <i>,  which  is  now  isotropic,  making  the  distribution  of 
0'  in  Eq . III.  3 immaterial.  The  final  result  is 


W(R)  = rJ^ P 


(pR)  F (p)  dp 


(m.6) 


This  equation  is  stronger  than  the  usual  relation  between  the  Hankel 
transform  and  two-dimensional  Fourier  transform  of  an  isotropic  function .[  104] . 
For  one  thing,  the  condition  of  isotropy  has  been  shown  to  be  unnecessary 
here;  furthermore,  when  the  characteristic  function  is  isotropic,  it  is  not 
only  a function  of  R alone,  which  is  true  in  general,  but  also  is  a function 
with  physical  meaning  as  a one -dimensional  Fourier  transform  (Eq . III.  4). 

The  extension  of  these  results  to  second-order  distributions  is  easily 
accomplished.  The  moments  are 


p (t)  = <x7y”  xP  y5  > = <R?‘'*'"r5'‘’^  cos"^  0 , sin"  0 , 

'^mnpq  ^1122  / av'l  2 1 1 


xcos^  0-  sir*^ 

2 2 /c 


av 


(ni.7) 

The  characteristic  function  can  be  put  in  the  form  corresponding  toEq.  (III.  3), 
namely , 


Pi  P2  ^2’  P2  ^2’ 


oo 


/•  .P 

(^  Pi ) (i  Pp)  w p w p 

- <‘^1  ®-2  (III. 8) 


M+P=0 


If  0j^  is  uniformly  distributed,  the  average  can  be  carried  out,  giving 


F^x  Y^Pl  dj,  pj  sind^;  P2COS  d2*  p2sind2;t)  = F2^(  p^ , P2‘.  t,  cos(d2~d  ^))  , 

(III. 9) 

that  is,  wherever  cos  oj^t  appeared  in  the  characteristic  function  for  V, 
cos  (d2~dj^)  will  now  appear. 


If  the  distribution  of  0^  is  not  rc.quired,  then,  as  before,  integrations 
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over  first  0j,  then  gives  the  analogue  to  Eq.  (III.  6) 

oo  oo  2u 


W2(Ri.R2,e.t)=BjR,J'  J J 


o o o 


d Pjdp^PiP^ 


1/2, 


■'o  (fp  1=^1  2*^2*  ) 


X (p  J,p  cos  , 


(in.  10) 


where 


® ~ ®2  " ®1  ' ^ “ ^2  ’ ^1 
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Appendix  IV. 

Clutter  Fluctuations  due  to  Movement  of  the  Scatter ers 

If  the  scatterers  in  one  resolving  volume  move  between  pulses,  e.g., 
are  blown  by  the  wind,  then  the  position  of  the  kth  scatterer  will  change 
in  the  time  t by  an  amount  d^^  = component  of  the 

scatterer' s velocity  along  the  radius  vector  to  the  radar;  t is  assumed 
short  enough  so  that  v^^  is  constant-  The  change  in  distance  will  alter 
the  time  of  arrival  of  the  pulse  reflected  from  the  scatterer  by  an  amount 


^k  = ^2k 


t'. 


Ik 

K 


2v,  t/c.  The  clutter  return  at  time  t is  then 
k z 


t , -t', , 2v,  t 

IK  K 


2v,_t 


cosja>^(t^+t)-o>,j.-^-  4^2k 
k=l  L 


(IV.  1) 


Because  the  start  of  each  transmitter  pulse  serves  as  a new  origin 
of  time,  the  time  difference  ~ appears  in  the  envelope  factor  only 

in  t^i  the  time  of  arrival  of  the  new  transmitter  pulse;  furthermore,  the 
phase  change  in  the  rapidly  varying  part  depends  on  the  transmitter  fre- 
quency, w-j,,  not  on  the  i-f  frequency,  From  this  equation  and  £q.  (3.  1) 

the  semi -invariants  are  found  to  be 


^CD 


\ (t) 

mn 


/ m n 
1*2 


h^(x)h"(x 


2vt . m 

) cos 

cf- 


(u)^7-x  - 


-00 


m , _ , t-VL  , ^ , V j V 

X cos  (<i)^Tx  + u,j.  — ^ + w^t  - ijj^)  dx^ 
where  x = t/  <f . 


2vt 


av 


(IV.  2) 


The  effect  of  moving  scatterers  may  be  illustrated  by  studying 
the  correlation  function  . We  have 


\ll(t) 

r(t)  

20 


= r (t)  cos  (w  t + a)  , 
o o 
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2 ioo^vt 


r„(t)  = r^(t)|  e 


i^\, 


>v, 


i|j  av 


- 4>)  > 

cos(u)^^  - ^)y 
r^(t)  = <a^a2>/<a^> 


tan  a = 

Here  i|<  = ■ 4*1  • 

and 


av 


av 


(IV.  3) 


h(x) h(x+ 


^oo 

i)dx/y  1 


-CD 


h (x)  dx 


CO 


This  equation  corresponds  to  the  one  first  found  by  Siegert,[  105]  except 
that  here  the  pulsed  nature  of  the  clutter  is  included  in  the  analysis,  ap- 
pearing through  the  correlation  function  of  the  pulse  envelope 

r j (2vt/c-r)  • 

In  particular  situations,  the  correlation  function  is  usually  simpler 
than  Eq . (IV.  3)  seems  to  indicate.  Almost  always  the  motion  of  the 
scatterers  varies  the  envelope  dependence  so  much  more  slowly  than  the 
phase  factor  that  the  envelope  correlation  function,  rj^(t),  may  be  taken  as 
unity  during  the  time  of  significant  variations  in  r^(t).  The  pulse  envelope 
correlation  function  will  drop  to  near  zero  when  the  argument,  2vt/cT’> 
is  near  unity,  i.e.,  a time  difference  comparable  to  the  pulse  duration. 

One  can  argue  that  up  to  2vt/c»j'ru  0-  li  this  correlation  function  may  be 
taken  as  one.  Then  the  velocity  average  becomes 


,2vt. 


2iu>>pvt 


W(v)dv  = 


-GO 


e ^ W(v)dv 


(IV.  4) 


With  a 1- jisec  pulse,  these  limits  are 


5.00 


, (velocities  in  m/sec),  so  that 


if  the  probability  of  velocities  exceeding  this  limit  is  sufficiently  small, 
the  limits  on  the  right-hand  side  can  be  made  infinite,  giving 


h ' c'T 


(^)  e 


ZitOrp  V t 


> 


2w_t 

= F (— ^-) 

V V C 


(IV.  5) 
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where  is  the  characteristic  function  of  the  velocity  distribution.  A 
reasonable  upper  limit  on  significant  velocity  values  is  v + 3 o^,  which 
for  "window"  is  about  2-5  m/sec.  Thus  up  to  times  of  the  order  of  one 
second,  the  correlation  function  assumes  the  simple  form  of  Eq . (IV.  5). 
The  argument  of  the  characteristic  function  is  then  sufficiently  large,  so 
that  the  correlation  function  r^(t)  is  negligibly  small. 

As  a specific  illustration,  consider  the  case  of  square  pulses  reflec - 
ted  from  clutter  with  a.  gaussian  velocity  distribution.  Then 


2iw.j,vt 


Ziw^vt 


<^h 


, 2v  t , 


> I H«‘- 


2;vt 
c r 


) e 


> 


(1  + - F 


2u)_  t 

<^'l 


4(r  cj  t 

1 2vt  . V T 
I 1 K 

c T 


.,222 
20*  w t 

V 1 


1 - vt  {6.7x10"^)  -i  o*^t^  (2.5x10"^) 


2 2 


800.0 


(IV.  6) 


Clearly,  the  exponential  will  ensure  the  effective  vanishing  of  the  whole 
expression  before  the  effect  of  the  envelope  correlation  function  can  be 
felt  at  all,  even  for  the  case  of  an  airborne  radar  measuring  ground  clutter 
(v/n,-500  mph  = 250  m/sec.  O'  /-^Im/sec). 

The  normalized  correlation  function  of  Eq.  (IV . 3)  is  thus 


I 

r (t)  = r (t)  F (-^)  F (-1,  1;  t) 
o a I V c 


(IV.  7) 


Furthermore,  the  change  in  amplitude  of  the  return  due  to  the  change 

in  range  is  so  small  as  to  be  negligible,  so  that  the  amplitude  correlation 

function  r (t)  depends  only  on  changes  in  aspect  of  the  scatterers  with  time. 
& 
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Where  the  aspect  of  the  scatterer  either  changes  very  slowly  or  does  not 
affect  the  anaplitude  of  the  return  (e.  g.  , "window"  sea  return),  then  r (t) 
is  also  unity.  The  phase  of  the  scattercrs,  is  usually  the  same  from 
pulse  to  pulse  so  that  the  correlation  function  reduces  to  the  magnitude  of 
the  velocity  characteristic  function  alone. 

Not  very  much  is  known  about  the  velocity  distribution  of  scatterer s. 
T.S.  George  [106]  has  considered  the  problem  of  a moving  radar  and  fixed 
clutter,  obtaining  fluctuation  values  through  numerical  integration  without 
explicitly  calculating  the  correlation  function.  Both  Hilst  [107]  and 
Fleisher  [108]  have  studied  the  inference  of  velocity  distributions  from 
experimentally  observed  correlation  functions;  however,  as  yet  the  time 
dependence  of  the  correlation  function  is  not  a settled  question.  In  par- 
ticular, all  of  the  three  types  of  time  dependence  considered  in  section  IV 
of  this  report  are  valid  possibilities  on  the  present  experimental  evidence. 
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Appendix  V 
Output  Power 

To  calculate  the  output  power  requires  fir>ding  a general  moment 
of  the  input  envelope,  as  Eq.  (4.4)  shows.  Using  the  nearly  gaussian 
envelope  distribution,  Eq . (2.26),  we  have 


2v+l  A R ■ 
— 


R^+A  ^ 
o 


The  integral  is  a generalization  of  Weber's  exponential  integral "[  1 09] 
well  known  to  give 


8i|j 


(V.2) 


where  ^( -v;  1; -p)  is  the  confluent  hypergeometric  function  and  p =A^/2i|i 
is  the  input  signal -to -noise  power  ratio.  The  indicated  differentiation  can 
be  carried  out  most  easily  by  using  the  series  for  the  hypergeometric 
function,  (Eq.  VI.  8),  to  give  the  final  result 

=r'(2+l)(24>)''  [ ^F^(-v;l;-p)  + — (v- 1) ^F ^(-v+2;  1 ;-p) +0(y ‘ ^)] 

3 ! 4* 

(V.3) 

For  the  d-c  power,  the  square  of  the  vth-moment  is  required,  which 
is  found  by  replacing  2v  by  v in  Eq.  (V.  3).  The  2vth-moment  is  slightly 
more  convenient  to  use  in  discussing  the  nearly  gaussian  correction  term 
so  the  next  paragraph  treats  only  the  total  power. 


The  nearly  gcvussian  deviation  in  total  power,  as  a function  of  v and 
p,  depends  on  v(v-l)  ^F  ^( -v  + 2;i  ;-p) . The  general  behavior  of  the  hyper- 
geometric function  can  be  found  by  examining  the  cases  of  integral  v, 
and  interpolating,  which  is  possible  because  the  hypergeometric  function 
is  a continuous  function  of  its  parameters.  A graph  of  these  functions  is 
shown  in  Fig.  17.  When  v is  greater  than  one,  both  v(v-l)  and  the  hyper- 
geometric function  are  positive,  so  the  correction  term  is  positive.  When 
V is  less  than  one,  v(v-l)  is  negative;  and  the  hypergeometric  function 
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changes  sign  at  a value  of  p greater  than  one.  Thus  the  correction  term 
is  negative  if  p is  less  than  one,  and  is  first  positive  and  then  negative 
if  p is  greater  than  one,  the  zero  of  the  correction  term  moving  toward 
V = I as  p increases. 
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Appendix  VI 

Evaluation  of  the  Output  Correlation 
Function 


The  total  output  correlation  function  is  given  by  Eq.  (4.  8), 


£(iz^) 


F,  (z,  ,z-)  . 

2 ' 1 2 sig 


(VI.  1) 


All  the  functions  appearing  here  are  known  except  for  the  characteristic 
function  of  the  sine  wave  signal,  which  can  be  evaluated  from  its  defini- 
tion as  an  average,  as  done  below.  To  form  an  ensemble  from  the 
determinate  sine  wave,  we  take  a random  phase  as  the  ensemble  par- 
ameter . 

X(t;d)  * A cos  (u)  t d)  , (VI.  2) 

o o 


where 


w(d)  * ^ «u:2« 

0 


Then 


F,(z,,z_)  . 

2 1 2 sig  2u 


’'2ir  . 


iz.A  cos  (w  t +^)  + iz_A  cos  (w  [t  ft]  + ^) 
lo  oo  2o  o'-  O'* 


d«( 


= J (A  k[z  , +z,  +2z,z,cosw  t) 
0 0*1  2 12  o 


oo 


c (-1)  J (A  z )J  (A  z )cos  nw  t , 
n n o 1 n o 2 o 


(VI.  3) 


n=0 


where  we  have  used  well  known  formulas  of  the  theory  of  Bessel  functions 
The  characteristic  function  of  the  noise  is 


^2^^1’^2^  = [1  + L]exp 


1 1 2 2 

+2i|ir^z,z,cosu)^t] 


(VI.  4) 


'2"2  o“l“2'  o 

where  L is  the  differential  operator  in  that  appear  s in  Eq . ( 2 
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The  correlation  function  can  be  separated  into  its  zonal  parts  by 
expanding  the  trigonometric  term  in  the  exponent  of  Eq.  (VI.  4)  as  a 
Fourier  series  and  combining  with  the  Fourier  series  of  Eq.  (VI.  3). 
The  result  is 


oo 


n=0 


oo 

€ COS  m w t 
m o 

m=0 


r(v+l) 

~ vV+1 

(iZj) 


^1  2 
e 


J (A  zj 
n o 1 


^2  2 
~Z~  ^2 

^ ^ JAAz) 

n o £, 


r,  , .m+n 
X[(-1)  € 


m+n 


I (4ir  z.z  ) + (-l)l”^“’^l  ,I  ,(\\iT  z.z,)]. 
m+n  ^ o 1 2 ' |m-n|  ]m-n|  ^ o 1 2 •* 

(VI. 5) 


This  equation  may  be  further  simplified  by  noting  that,  for  the 
receiver  problem,  only  the  low  frequency  zone,  m=0,  is  required.  Then 
the  two  contour  integrals  may  be  separated  by  expanding  the  modified 
Bessel  function  in  a power  series,  giving 


R(t)  = [l  + L] 


,'^^o.n+2k 

-irnuV-ri M i,(A  .zJH  .(A  ,z,) 

A k!(k+ny!  n,k  o 1 n,k  o 2 


n=0  k=0 


(VI.  6) 


4'i=^1>2=4^ 


where 


_pr  (v+1) 


vb.z . 

-iJ- 


H AA  ,z.)=e.'.  / dz.  e (A  z.)  ^ 


and  the  low-frequency  correlation  function  R(t)  appears.  The  H , func- 

n j K 

tions  may  be  evaluated  by  taking  as  the  contour  the  real  axis  indented 
downward  at  the  origin  (see  sec.  4 and  App.  Ill  of  ref.  100). 


H .(A  ,z.)=eQii^[i 

n,k  o j 


n+2k-v 

X 


1 

e 


J (A  x)dx 
n o 


provided  a + k^  v/2,  so  that  the  integral  around  the  indentation  vanishes. 
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The  real  integral  is  Weber's  first  exponential  integral  [ 110]  , so  that  the 
final  result  is  n v n+2k 

.T  X , ,,n+k  ^ 

H ^(A  ,z  ) = (-1)  -=T F (-.y+n+k;n+l;-p  ) 

^ (VI. 7) 

2 

In  this  equation,  p^  = 1^1  confluent  hypergeometric  func- 

tion. The  generalized  hypergeomdric  function  is  defined  formally  as 


F (a . , « • . a J b . , . . . b ^x^= 
p q 1 pi  q 


oo  (a- ) 

— 7 . , ^ n n 

1=  1 


X 

n ! 


(VI.  8) 


n=  o J7  (b . ) 
i=l  " 


where 


I 1 = 0 


In  this  report,  only  the  three  functions  jFj^,  2^0  needed. 

Equation  (VI.  7)  iJ  an  integral  function  of  n + k - v/2  and,  therefore,  by 
analytic  continuation  k^‘^o‘^j^  equals  the  right  side  of  Eq.  (VI.  7)  with- 
out any  restriction  on  v.  The  coi  relation  function  is  thus 
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oo  OO 


‘n'PlPz’ 


n 

7. 


L ^ V-^kinl^k-tn)!  \ 

n=o  k=o  ' ^ 


! \ 


n-tZk 

2 2v  ~T~ 


' Tn+ic 


X iFi(-^-tn-tk;n-t  1;  -Pf)  iFi(-.^-tn-tk;n  + 1;  -p^) 


(VI.  9) 


47j-iji2=4^ 


By  expanding  the  hypergeometric  functions  and  summing  over  k, 
an  alternate  form  may  be  found,  viz.. 
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R(t)  = [1  + L] 


p V^v+1) 
4r^(7+l) 


V 
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. , n . . n n 

i+T  j + 

oo  cr>  oo  , , , .1+1  2 

‘^(-1)  -Pi  P2  r 
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1 ✓ , 2 2x7 
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+ 2 

di  r 

+ -^+j+n;n+  1; 


(VI.  10) 


When  p is  small,  this  last  is  the  better  form,  since  it  is  a power 
series  in  p.  For  large  p,  the  ordinary  asymptotic  expansion  of  the 
confluent  hypergeometric  function  may  be  substituted  into  Eq.(VI.S)  to 
give  a formal  asymptotic  expansion  of  the  correlation  function 

^F^(a;c;-x)c^P^^^j^^  x"^  ^^^(a,  1+a-c;^)  + O (e"^)  , (VI.  11) 


|xl>>  a , jx|»  c 


R(t)c^r[l  + L] 


pV(v+l) 

4r^(Jfl) 
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oo  CD 
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t -J  ->  r n+2k 

n , Vv2  I V)2  _o. 

k! (k+n) ! 2'k  2'n+k  ' p ' 


n=o  k-o 


^ 2^0^'  2"'’'^''’^’  ■ 2 "’’^'p,  ^ 2^0^"  2'*’ " 2'*’ ^ 

X w 


Notice  that  as  p tends  to  infinity,  the  correlation  function  reduces  to 
that  for  the  signal  alone,  as  it  should;  however,  the  asymptotic  nature  of 
the  equation  is  not  assured  ac  present,  since  Eq.  (VI.  11)  is  not  valid  for 
values  of  n and/or  k comparable  to  or  greater  than  p.  In  these  cases,  by 
using  the  asymptotic  expansions  appropriate  to  n>>  p,  k;  k^>^  p,  n;  n,k]^^  p 
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n = p;  k = p;  n = k = p;  [ 1 1 1]  it  can  be  shown  that  the  leading  term  in  the 
asymptotic  expansion  is  of  exponential  order.  Thus  the  terms  which  are 
neglected  in  using  a finite  part  of  Eq.  (VI.  12)  are  actually  of  even  smaller 
order  than  the  equation  indicates  so  that  the  asymptotic  expansion  is  valid, 
with  the  restriction  that  p must  be  larger  than  the  largest  value  of  n or  k 
appearing  in  the  finite  approximation.  In  the  event  that  v is  an  even  integer, 
all  the  series  involved  terminate,  and  the  left-and  right-hand  sides  ofEq.(VI.  12) 
are  equal. 

The  remaining  step  is  to  evaluate  the  derivatives  contained  in  the  operator 
L.  The  strong  signal  case,  Eq.  (VI.  12),  contains  only  positive  powers  of  the 
quantities  to  be  differentiated,  so  that  the  result  is  simply  obtained  using  the 
formula 


3—  .,F.(a,b;cx)  = abc  TF.(a  + 1 , b + 1 ; cx) 
ox  c V U 


(VI.  13) 


The  result  is.  in  abbreviated  form. 
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(VI.  14) 


where 


2^n+k+a^“  2^k+a  2*  0^  ‘ ‘Z"*' ^ ‘ 'Z'*' ^ p ^ 

The  weak  signal  case  is  more  complicated,  but  with  the  analogue 
of  Eq.  (VI.  13)  for  the  function,  and  the  additional  formula 

^ [x^  2^j(a,b;c;x)]  = ax^“  ^ ^(a+ 1,  b;  c;  x)  , 


the  result  corresponding  to  Eq.  (VI.  14)  is,  also  in  abbreviated  form. 
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- oo  CD  oo  £ , , ,m+n  k k+m+n 
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where 


3(a,b,c)  = (-l)k+„+a<-l'k+m+b  j(  - |tk+„+a.-|+k+m+b;  k+Uc:  . 


When  V is  an  integer,  the  hypergeometric  functions  reduce  to  polynomials 
if  V is  even,  or  to  complete  elliptic  integrals  if  v is  odd. 
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Appendix  VII 

Effect  of  Nearly  Gaussian  Statistics  on  the  Output 
Correlation  Function 

To  determine  whether  the  correlation  in  the  output  is  increased  or 

decreased  by  nongaussian  noise  as  compared  to  gaussian  noise,  we  study 

various  special  cases  of  r ^(t). 

^ out 

When  p is  large,  we  may  use  the  strong  signal  expression  for  the 
correlation  function,  Eq.  (4.  13),  to  obtain  finally 


AjiW+AjjU) 


r 1 r - +A,3(t 


- 


(vn.  1) 


For  the  exponential  and  gaussian  models,  we  see  from  Table  4.  1 that  the 
semi -invariant  expressions  in  the  inner  braces  are  negative  so  that  when 
V is  greater  than  two,  the  output  correlation  is  reduced.  When  v is  quite 
small,  the  second  term  alone  is  important  and  again  the  output  correlation 
is  reduced.  Near  v = 1,  the  sign  of  the  change  in  correlation  varies  be- 
cause the  output  correlation  functions  cross  each  other.  The  boundaries 
between  this  region  and  those  in  which  the  nongaussian  output  correlation 
is  less  lie  in  the  ranges  v = 0.4-.6,  l.Z-1.4for  both  models . 

The  boundary  is  not  sharply  defined  since  inside  the  region  the  gaussian 
and  nearly  gaussian  curves  cross  each  other. 

For  the  linear  model  the  nongaussian  part  of  Eq.  (VII.  1)  vanishes. 

We  must  then  star;  again  from  the  general  strong  sigiml  result,  Eq.  (VI.14), 
_3 

including  p terms  specialized  for  the  linear  model.  Carrying  this 


procedure  through  gives 


r (t)  = r Jt) 
out  out  gauss 


[7v^  - 54v^+ lZ6v^  - 1 28v  + 56]  + 0(p  ^,-y  ^ 


(VII.  Z) 
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The  difference  between  gaussian  and  nongaussian  noise  is  of  smaller 
order  here  than  for  the  other  two  noise  models.  The  polynomial  has 
two  real  roots,  when  v equals  1,55  and  4.54,  between  which  points  the 
nongaussian  result  has  more  correlation. 

When  no  signal  is  present,  Eq.  (VI.  15),  which  contains  hyper - 

geometric  functions,  must  be  used.  The  behavior  of  the  correlation 

function  as  a function  of  v is  then  still  too  difficult  to  predict  for  all 

times  L;  however,  by  expa.nding  the  hypergeometric  functions  in  powers 
2 2 

of  r^  (t)  and  1 - r^  (t),  the  change  in  correlation  can  be  found  separately 
for  large  times  and  small  times.  Since  the  two  results  agree,  we  have 
confidence  iu  the  indicated  answers. 


For  large  times,  the  input  correlation  function  is  small  and  the; 
expansion  of  the  hyper  geometric  function  around  the  origin  is  appropriate 
(Eq.  (VI,  8)).  The  normalized  output  correlation  functions  for  the  noise 
models  are  then 


, ^2  ■ ^ „ r -2 

= r^  I 1 - — + O [y  , 


out' 


2 , 2 


-1 


r -4.  2 ^ ! a'\ 

[y  ,v  r^  /4] 


-2  2 2, 


, exponential  model  (VII.  3) 


^^^(t)  = r^  I 1 — ^ (3-2/^^r^) + 0[y  ,v  r^  /4]>  gaussianmodel  (VII.  4) 


r ,(t)  = r 1 + 
out  O I 


+ 0[y'^.  v^/4] 


linear  model  (VII.  5) 


These  equations  indicate  that  nongaussian  noise  of  the  linear  type  is  al- 
ways more  correlated  than  the  corresponding  gaussian  noise,  with  identi- 
cal results  for  the  exponential  and  gaussian  models  when  v is  less  than 
two,  and  opposite  results  for  v greater  than  two. 


For  small  times,  the  expansion  of  the  hyper  geometric  function  in 
2 

powers  of  1 - r is  needed.  [ 1 12]  . 
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F(a,b;c;z) 


_T(c)r(c-a-b) 

p(c-a)f'(c-b) 


F(a,b;a+b-c+  1;1  -z) 


F(c]f^(a+b-c) 
p(a)p  (b) 


(1-z)^  ^ ^F(c-a,c-b;c-a-b+l;l-z) 


(VII.  6) 


Integral  values  of  c-a-b,  which  correspond  to  integral  values  of  v, 
are  exceptional;  however,  the  modifications  of  Eq.  (VII.  6)  that  are  then 
needed  when  v is  greater  than  one  are  of  higher  order  than  the  terms  ap- 
pearing in  Eqs  . ( VII.  7-9)[  1 13]  . For  the  linear  detector  the  results  are 
already  known  (Table  4.  3)  so  that  a special  treatment  cf  this  case  here  is 
unnecessary. 


The  normalized  output  correlation  functions  are 
exponential  model : 


r x(t)  r ,.(f)  ^ 1 - 

cut  out  gauss 


% 5-  -^-[30-4-.-  . +2v-4)] 

p^(^+i)-|2  riv+i) 


2,v 


r^(j+ 1) 


34,' 


.r 

[‘-fTTF- 


gaussian  model; 


r (t)  = r (t) 
out  out  gauss 


linear  model: 


-2 


2v2, 


r (v+1) 


+ 0(Y  “(1-r  ")“) 


(VII.  7) 


2,,  2 


[1  . 

-rr^J 


r(T+ir(T+2)' 


T (v  + l) 


+ 0(y^,{\-rl)^'p 


(vn.8) 


r ,.(t)  = r (t) 
cut  out  gauss 


h40U-r‘) 

L'twtJ 


V 

~ ^ 
V - . 


r^/v  + i)  1 

-,V  .,2^'  '2  ' / 3 3 2^_ 

^^2'^^  ■‘TTvTiy  +2v-2) 


-2  22, 

+ — a- +0(Y 


6i|i 


if  p(v+i) 

-p(v+l)J 


(VII.  9) 
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Appendix  VIII 

The  Effect  of  Nearly  Gaussian  Statintics  with  Finite  Averaging 

When  the  signal-to-noise  ratio  is  large,  the  output  correlation  func- 
tion is  comparatively  simple,  Eq.  (4.  13),  so  that  the  variance  for  general 
V can  be  found,  following  Eq.  (4.  15), 

exponential  model: 


“7 

X 


4pY  2'2 

gaussian  model; 


+ l)(e'^*+3x-  1]  + 0(p'^,-,'^)l 


2 

G X ^ p J 'VI 

k. 


pV2 


+ v{v-Z)  ^ ^ 0(p‘^,  y'^) 

pyV^ 


(vm.  1) 


2p'/Z’ 
2 

e-^  -1] 


SpyVtt’ 


^ -1] +o(p-^  y‘^) 


(VIII.  2) 
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linear  model: 


c^,|[[2  + 2<^'+j^[2(|)^v(v-2)  + (|-1)2]]  [3-x] 

2 

+ [(^)*^  + (^-  ^)^]  [3  - 2x  +^]  +0(p  ^ , x<  1; 


■'■  [4x-l]  +0(p'“,  y“^)  x^  1; 


where 


= pT 


i>{x)  = ^ 

vir 


e"^  dt 


K = 

4r^(4  + 1) 


When  the  signal  is  large,  the  nongaussian  terms  are  all  of  higher  order 

than  p ^ and  so  will  not  appear  in  I to  the  order  used  above.  Thus 

2 _2 

the  effect  of  the  nearly  gaussian  statistics  on  O'  /I  is  entirely  fixed  by 
the  terms  appearing  in  Eqs.  (VIII.  1-3).  For  large  x,  the  variances  re- 
duce to 

expoaential  model: 

2 2K  r.,  , r ^ y , V , ^ , V , 2 . V , . 2, 


=-;rp2p  [^2*2-  S-2'  ) 


, -1  - 2 -2, 

+ U(x  , p , y ) 
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gaussian  model: 


2 kK^ 


= - 
G X 


P’jfV^TT 

linear  model: 


^ r/'^\2  , 4 V -V  n , 1 /V  i\2-|  -1  -2  -2. 


(vm.5) 


+ + +0(x-\p-^.  y‘-) 


(vm.6) 


The  nongaussian  part  of  all  three  of  these  is  positive  for  large  and 
small  V and  negative  for  v near  one.  The  two  values  of  v at  which  the  non- 
gaussian part  of  the  variance  changes  sign  are  for  the  exponential  model 
0 . 33  and  1.42,  for  the  gaussian  model  0 . 38  and  1 . 33,  and  for  the  linear 
model,  0.46  and  1.26. 

For  the  no-signal  results,  it  is  most  convenient  to  start  from  an 
alternate  expression  for  the  relative  variance. 


^2 


R(t) -I  _ R(t)  - R(oo)  R(0)  - R(oo) 
■ R(0)  - R(cx3)  ’ STooT 


= r (t) 
out 


(vm.7) 


This  form  enables  us  to  use  the  results  already  found  for  the  behavior 
of  For  noise  alone,  from  Eqs.  (4.5,  4.6),  we  have 


^o-^d-c  r(vti) 

— p 


d-c  r^+i) 


1 - ..  + + 0('»— ) 


(VIII.  8) 


so  that  the  nongaussian  part  of  the  power  ratio  is  always  positive.  The 
relative  variance  is  then 


J /^0-Pd-r\2  f 
V “"d-C  /V 

' '-'O 


(VIU.9) 
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in  which  1 - t/*r  is  always  positive. 


Since  is  always  greater  for  nearly  gaussian  noise  for  the 

linear  model,  it  follows  immediately  that  the  nongaussian  variance  is 
always  greater.  When  v is  less  than  tv/o,  the  same  is  true  of  the  ex- 
ponential and  gaussian  models;  however,  for  these,  when  v is  greater 
than  two,  is  less  for  nongaussian  noisc,  so  the  answer  is  not  im- 

mediately obvious.  For  the  two  models  with  r^^  small,  from  Eqs .(VII.  3,4) 

exponential  model: 


R(t)-r^_  2 P(V+1)  f,  ^ I [ 


(Vm.  10) 


gaussian  model; 


R(t)  • P_  j.  2 r(v-n)  I 


1 - 


p(^+l). 


' '7^^'  12  r ! 


+ 0(y‘^,  r^^) 


(vm.  11) 


The  nongaussian  part  of  both  of  these  are  now  positive  for  all  values 
2 

of  V.  When  i - r^  is  small,  from  Eqs.  (VII.  7,8),  we  have 
exponential  model: 


,ti  J^U 

■jZ  out'  gauss  ^2^1/4.!) 

2 


V 

V - 


32t{l-[r^(^)+l)/hv+l)]} 

[•'(4-3  [1  -I-/])  + 4 r^^  + E^^^|^4r^^  - 2v(  1+r^^)  + v ^ 1 - r ^^)j 


+ 0(Y'^(l-r^^)^) 


(Vin.  12) 
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gaussian  model; 

R(t) -r^-^  [n  

out  gauss  32yV;?11  -r^(^+l)/r(v+l)] 


f 21  / r^T+ih 

x[(4-2[l-rj)(l 


+ o(y‘^.(i- 


2.2v 
^o  ) ) 


(vm.  13) 


These  are  also  larger  when  the  noise  is  nongaussian,  so  that  finally  we 
have  the  result  that  the  variance  is  always  larger  with  nongaussian  noise 
than  with  gaussian  noise  when  no  signal  is  present. 


This  same  approach  can  be  used  to  find  the  change  in  the  variance 
when  nongaussian  noise  is  present  and  the  signal  is  strong.  The  power 
ratio  is  now 


o d-c 

“PT 

d-c 


2 

V 


. ^2v^+3v+8  A 


40  7v_  - 12v+ 14^Q^p-2^^-2^  1 ^ (vill.  14) 


6pi}< 


The  nongaussian  part  is  negative  between  v = 0.45  and  v = 1.26.  Since  the 

-2 

nongaussian  part  of  linear  noise  model  is  of  order  p , the 

power  ratio  alone  will  determine  the  sign  of  the  change  in  the  variance. 

So  with  these  two  noise  models,  the  zeros  of  Eq  (VIII.  14)  correspond 
rather  closely  to  those  of  so  that  when  vis  small,  the  nongaussian 

variance  is  more;  when  v is  near  one,  it  is  less;  and  when  v is  two  or 
more,  it  is  greater  again.  Finding  more  exact  boundaries  by  this  method 
is  not  worth  while  since  they  are  already  known  from  the  first  part  of  this 
appendix. 
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